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ABSTRACT 



Radiative corrections to processes that involve the production and subsequent decay of un- 
stable particles are complex due to various theoretical and practical problems. The so-called 
double-pole approximation offers a way out of these problems. This method is applied to the 
reaction e~^e~ — > W^W~ —>■ 4 fermions, which allows us to address all the key issues of dealing 
with unstable particles, like gauge invariance, interactions between different stages of the re- 
action, and overlapping resonances. Within the double-pole approximation the complete 0{a) 
electroweak corrections are evaluated for this off-shell ly-pair production process. Examples of 
the effect of these corrections on a number of distributions are presented. These comprise mass 
and angular distributions as well as the photon-energy spectrum. 
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1 Introduction 



In order to test the Standard Model (SM) of electroweak interactions, millions of Z bosons 
have been produced and studied at LEPl High-precision measurements of the Z-boson 
parameters have been performed by measuring 

e+e- ^Z^ff (1) 

and comparing the results with the theoretical predictions. This requires theoretical precision 
calculations for the reaction 

e+e- ^ //, (2) 
involving all lowest-order diagrams and the radiative corrections (RC) to them . 

Measuring the parameters of the W boson provides further tests of the SM. In particular, the 
accuracy of the VT-boson mass {M\y) has to be improved, since it can give indirect information 
on the Higgs sector and on physics beyond the SM. At LEP2, W bosons can be studied in the 
I^-pair production reaction 

e~^e^ W^W~ 4 fermions. (3) 

Besides the mass Mw also the Yang-Mills form of the triple gauge-boson couplings (TGC) 
ZW^W~ and '-)W^W~ can be tested in reaction (|^), where they directly manifest themselves 
in the lowest-order cross-section. Deviations from the SM Yang-Mills couplings can be searched 
for most effectively by a detailed investigation of angular distributions [Q. 

Just like in the case of the Z boson, the precise determination of the IV-boson parameters 
requires both accurate experiments and accurate theoretical predictions. Since the statistics at 
LEP2 is much smaller than at LEPl, the theoretical calculations do not have to be as precise 
as those for LEPl. However, just the Born prediction for process (^, involving only three 
diagrams, is not sufficient 0. Therefore, one would like to have an idea how the cross-sections 
for the four-fermion process 

e^e~ —>■ 4 fermions (4) 

are affected by the inclusion of the remaining lowest-order diagrams and the RC to the complete 
set of lowest-order diagrams. 

The question of the complete lowest-order calculation of process @ has been studied in the 
literature |[. Roughly speaking, those diagrams that contain a single W boson are a factor 
of 0(Tw/Mw) smaller than those for the VT-pair production process (|^). Diagrams that do 
not contain a W boson at all are down by 0(r^/M^). Note, however, that some diagrams 
that do not contain two W bosons can nevertheless be large, e.g. as a result of the exchange 
of almost real photons. Besides such special cases, the lowest-order "background" diagrams in 
(^), i.e. the non-M^-pair diagrams, will give at most a correction of 0{rw/Mw) to the Born 
cross-section of process @. 

The RC to process (|]) are a priori of 0{a), i.e. they are genericly of the same order as 
the lowest-order background diagrams. In analogy to the lowest-order case, the RC to the 
background diagrams are at most of 0{aT]y/M]y). 
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In view of the above estimates, the most relevant corrections to the Born cross-section of 
process are divided into two classes: 

1. 0{a) RC to the VT-pair process (|). 

2. 0(Tw/Mw) lowest-order background contributions from the full four-fermion process (H). 

At present the first class of corrections has not been fully studied in the literature. What has 
been discussed quantitatively and qualitatively is merely a subset of 0{a) effects: 



1.1. Initial-state radiation (ISR). Since ISR corrections are enhanced by collinear-photon ef- 
fects, they are large and even higher-order contributions should be taken into account 

i-i- 

1.2. Final-state radiation (FSR). Similarly FSR can be sizeable, in particular for distributions 
in the fermion-pair invariant masses. This QED effect has recently been pointed out in 
the literature [0. Again higher-order corrections are non- negligible for these W line-shape 
distributions. 

1.3. Whereas the effects 1.1 and 1.2 are enhanced by logarithms originating from collinear 
photons, also non-enhanced 0{a) effects have been studied. One of those belongs to the 



class of the so-called non-factorizable corrections P]-Il0[, i-e. corrections that at first sight 
do not seem to have two overall W propagators as factors. As such one would expect 
these corrections to be suppressed by additional powers of Tw/Mw- In the special case 
of semi-soft photonic corrections, however, the suppression does not take place and an 
0{a) correction survives. This 0{a) correction has been calculated and turns out to be 
relatively small. In the vicinity of the VT-pair threshold also the Coulomb interaction 



between the unstable W bosons has been studied in great detail |]TT| . The corrections are 
relatively large, but higher orders are not required. 

1.4. The effects 1.1-1.3 are all of QED origin. Complete electroweak RC to process (|]) have 
not been applied yet, but some attempts have been made to take into account the dom- 
inant effects. An overall effect of electroweak corrections has been considered by using 

as coupling constant instead of a. From other calculations it is known that a 
parametrization of the lowest-order term reduces the size of the one-loop non-photonic 
RC considerably. In addition to this overall effect, the full electroweak RC to stable 
ly-pair production [jl2|-[jl4| and on-shell W decay [T^ are already known for quite some 



time. 



The purpose of this paper is twofold. In the first place we will present a complete quan- 
titative evaluation of 0{a) electroweak RC to ly-pair-mediated four-fermion production. To 
this end we calculate all 0{a) factorizable corrections and add them to our previous results on 
non-factorizable corrections. In this way the gap in the above list of corrections 1.1-1.4 will be 
filled, and the exact size of all 0{a) corrections to reaction will be known. In the light of 
the physics motivation given above, this practical result is clearly wanted. 
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Secondly, this paper discusses in detail the so-called pole scheme. This scheme offers a way 
to avoid theoretical problems associated with the gauge-invariant treatment of reactions that 
involve the production and subsequent decay of unstable particles. The study of the W^-pair 
case serves as a relevant example of this method and can be used as a guideline for other cases. 
Since there are more unstable-particle production processes to be studied at future accelerators, 
the relevance of the presented study goes beyond the TV-pair case. 

The outline of the paper is as follows. In Sect. ^ various gauge-invariance issues will be 
discussed. In Sect. ^ the pole-scheme treatment is described, with special emphasis on the 
so-called double-pole approximation (DPA), since it is the basic ingredient for our calculation 
of the 0{a) RC. The comparison with the exact lowest-order evaluation in Sect. ^ gives an 
estimate of the accuracy of this approximation. The discussion of the RC in the DPA for the 
IV-pair production process (|^) is presented in Sect. |^, whereas the corresponding numerical 
results can be found in Sect. In Sect. ^ we give the conclusions of our study. 



2 Gauge-invariant treatment of unstable gauge bosons 
2.1 Lowest order 

The above-described processes, with or without RC, all involve fermions in the initial and 
final state and unstable gauge bosons as intermediate particles. Sometimes a photon is also 
present in the final state. If complete sets of graphs contributing to a given process are taken 
into account, the associated matrix elements are in principle gauge- invariant, i.e. they are 
independent of gauge fixing and respect Ward identities. This is, however, not guaranteed for 
incomplete sets of graphs like the ones corresponding to the off-shell PT^-pair production process 
(^). Indeed this process was found to violate the SU{2) Ward identities 0. 

In addition, the unstable gauge bosons that appear as intermediate particles can give rise 
to poles — M^) in physical observables if they are treated as stable particles. This can be 
cured by introducing the finite decay widths for these gauge bosons. In field theory, such widths 
arise naturally from the imaginary parts of higher-order diagrams describing the gauge-boson 
self-energies, resummed to all orders. However, in doing a Dyson summation of self-energy 
graphs, we are singling out only a very limited subset of all possible higher-order diagrams. It 
is therefore not surprising that one often ends up with a result that violates Ward identities 
and/or retains some gauge dependence resulting from incomplete higher-order contributions. 

Since the latter gauge breaking is caused by the finite decay width and is, as such, in 
principle suppressed by powers of F/M, one might think that it is of academic nature. For LEPl 
observables we indeed know that gauge breaking can be negligible for all practical purposes. 
However, the presence of small scales can amplify the gauge-breaking terms. This is for instance 
the case for almost real space-like photons [|^, [l^ or longitudinal gauge bosons (Vl) at high 
energies |]T3, involving scales of Olp^ /E"^) for B = 7, Vl- The former plays an important role 



in TGC studies in the reaction e'^e e Ueud, where the electron may emit a virtual photon 
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with an invariant mass as small as m^. The latter determines the high-energy behaviour of 
the generic reaction e^e~ — ^ 4 fermions. In these situations the external current coupled to the 
photon or to the longitudinal gauge boson becomes approximately proportional to . Sensible 
theoretical predictions, with a proper dependence on and a proper high-energy behaviour, 
are only possible if the amphtudes with external currents replaced by the corresponding gauge- 
boson momenta fulfil appropriate Ward identities. 

So, how should one go about including the finite decay widths? The simplest approach 
is the so-called "fixed-width scheme", involving the systematic replacement — My) —>■ 

~ My + iMyTv), where Fy denotes the physical width of the gauge boson V with mass 
My and momentum pv- Since in perturbation theory the propagator for space-like momenta 
does not develop an imaginary part, the introduction of a finite width also for py < has no 
physical motivation and in fact violates unitarity, i.e. the cutting equations. This can be cured 
by using a running width iMyTviPy) instead of the constant one iMyTy ("running-width 
scheme"). However, as in general the resonant diagrams are not gauge-invariant by themselves, 
the introduction of a constant or running width destroys gauge invariance. 

A truly gauge-invariant scheme evidently has to be a bit more sophisticated than this. It 
should be stressed, however, that any such scheme is arbitrary to a greater or lesser extent: 
since the Dyson summation must necessarily be taken to all orders of perturbation theory, and 
we are not able to compute the complete set of all Feynman diagrams to all orders, the various 
schemes differ even if they lead to formally gauge-invariant results. Bearing this in mind, we 
need besides gauge invariance some physical motivation for choosing a particular scheme. In 
this context two options can be mentioned. The first option is the so-called "pole scheme" |TP|. 



In this scheme one decomposes the complete amplitude by expanding around the poles. As the 
physically observable residues of the poles are gauge- invariant, gauge invariance is not broken 
if the finite width is taken into account in the pole terms oc l/(py — My). Note that the 
leading terms in such an expansion play a special role in view of their close relation to on-shell 
production and decay of the unstable particles. This point will be explained in more detail in 
the following sections. 

The second option is based on the philosophy of trying to determine and include the minimal 
set of Feynman diagrams that is necessary for compensating the gauge violation caused by the 
self-energy graphs. This is obviously the theoretically most satisfying solution, but it may cause 
an increase in the complexity of the matrix elements and consequently a slowing down of the 
numerical calculations. For the gauge bosons we are guided by the observation that the lowest- 
order decay widths are exclusively given by the imaginary parts of the fermion loops in the 
one-loop self-energies. It is therefore natural to perform a Dyson summation of these fermionic 
one-loop self-energies and to include the other possible one-particle-irreducible fermionic one- 
loop corrections ( "fermion-loop scheme" ) [0, |T^. For the LEP2 process e^e~ 4 fermions this 
amounts to adding the fermionic corrections to the triple gauge-boson vertex. The complete 
set of fermionic contributions forms a gauge-independent subset and obeys all Ward identities 
exactly, even with resummed propagators |jl8| . 



The above arguments, although general, apply in particular to tree diagrams. Therefore an 
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additional discussion for RC is required. 
2.2 Radiative corrections 

The next question that should be addressed involves the interplay between RC and gauge in- 
variance. After all, the RC are indispensable for coming up with reliable theoretical predictions 
for the LEP2 process e^e~ — ^ 4 fermions. 

As far as real-photon corrections are concerned, not much changes as compared to the 
lowest-order case. Still both the pole scheme and fermion-loop scheme yield gauge-invariant 
results. However, in the fermion-loop-scheme treatment of the process e^e~ — >■ 4/7 the full 
set of fermionic corrections to the quartic gauge-boson vertex emerges. This evidently is too 
much complexity for a tree-level calculation. The pole-scheme, with its close relation to on-shell 
subprocesses, remains relatively simple. As we shall see later on, some subtleties arise when 
photons are radiated from a virtual W boson, because this W boson may give rise to two poles. 

The implementation of the one-loop RC adds an additional level of complexity by the sheer 
number of contributions (10^-10'^) that have to be evaluated. By employing a gauge-invariant 
lowest-order finite-width scheme it is possible to cover the most important electroweak effects, 
like running couplings and leading QED corrections (see previous section), which are controlled 
by factorization theorems. However, there is still the question about the remaining corrections, 
which can be large, especially at high energies ^, pO[ . 

In order to include these corrections one might attempt to extend the fermion-loop scheme. 
At present this solution is not yet workable in view of the fact that a gauge-invariant inclusion 
of the one-loop corrections to the decay width in turn requires the inclusion of (the imaginary 
parts of) some two-loop corrections. Moreover, the number of one-loop contributions that have 
to be evaluated remains large. 

As a more appealing and economic strategy we discuss in the next section how the RC can 
be calculated in an approximated pole-scheme expansion. 

3 The pole scheme in double-pole approximation 

As mentioned above, the pole scheme consists in decomposing the complete amplitude by ex- 
panding around the poles of the unstable particles. The residues in this expansion are physically 
observable and therefore gauge-invariant. The pole-scheme expansion can be viewed as a gauge- 
invariant prescription for performing an expansion in powers of YyjMy. It should be noted 
that there is no unique definition of the residues. Their calculation involves a mapping of off- 
shell matrix elements with off-shell kinematics on on-resonance matrix elements with restricted 
kinematics. This mapping, however, is not unambiguously fixed. After all, it involves more 
than just the invariant masses of the unstable particles and one thus has to specify the variables 
that have to be kept fixed in the mapping. The resulting implementation dependence manifests 
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itself in differences of subleading nature, e.g. 0{Vy / My) suppressed deviations in the leading 
pole-scheme residue. In special regions of phase space, where the matrix elements vary rapidly, 
the implementation dependence can take noticeable proportions. This happens in particular 
near phase-space boundaries, hke thresholds. 

In order to make these statements a bit more transparent, we sketch the pole-scheme method 
for a single unstable particle. In this case the Dyson resummed lowest-order matrix element is 
given by 



Pv 



E 



„=o 




W{pI,uj) 



UJ 



1 



+ 



W{p 



1 



Pv 



Ml 



(5) 



where T,v{p 



Ly + EviP 

is the unrenormalized self-energy of the unstable particle V with momentum pv 
and unrenormalized mass My- The renormalized quantity is the pole in the complex py 
plane, whereas Z(M^) denotes the wave-function factor: 

- M^ + tviM^) = 0, Z{M^) = l + t'y{M^). (6) 



The first term in the last expression of Eq. represents the single-pole residue, which is closely 
related to on-shell production and decay of the unstable particle. The second term between 
the square brackets has no pole and can be expanded in powers of py — M^. The argument 
u denotes the dependence on the other variables, i.e. the implementation dependence. After 
all, the unstable particle is always accompanied by other particles in the production and decay 
stages. For instance, consider the LEPl reaction e~^e~ — *• //. In the mapping p| —>■ one 
can either keep t = [p^- — pf)"^ = — — cos0)/2 fixed or cos^^. In the former mapping 
cos6'poic is obtained from the on-shell relation cos^poic = 1 + 2t/M^, whereas in the latter 
mapping tpoic = — M^(l — cos d)/2. It may be that a particular mapping leads to an unphysical 
point in the 'on-shell' phase space. In the present example tpoie will always be physical when 
cos 6* is kept fixed in the mapping. However, since | cos^poiel > 1 for t < — ReM^, it is clear 
that mappings with fixed Mandelstam variables harbour the potential risk of producing such 
unphysical phase-space points.^ This can have repercussions on the convergence of the pole- 
scheme expansion. Therefore we choose in our calculations only implementations that are free 
of unphysical on-shell phase-space points. 

It should be noted that the mass and width of the W and Z bosons are usually defined in 
the so-called on-shell scheme: 



M^-M^ + Rety{M^) = 0, 



Zos{M^) = l + Ref:'y{M^), 



MyT 



yi y 



lmty{M^ 



v> 



■ (7) 



Both schemes can be related according to (see e.g. Ref. fTSlI ): 



{Ml - iMyV^ 



1 - 



+ 



M3 



^In the resonance region, \p\ — M'^\ <C the unphysical 'on-shell' phase-space points occur near the edge 

of the ofF-shell phase space, since t < — ReM^ requires cos 6* w — 1. 
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production 
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Figure 1: The generic structure of the factorizable VT-pair contributions. The shaded circles 
indicate the Breit-Wigner resonances, whereas the open circles denote the Green functions for 
the production and decay subprocesses up to 0{a) precision. 



As we are aiming for 0{a) precision in our study, the differences between both schemes can be 
ignored. For the same reason ipyTv/Mv can be replaced by iMyTy, since the difference only 
induces C(a^) corrections to the cross-sections. 

The at present only workable approach for evaluating the RC to resonance-pair-production 
processes, like ly-pair production, involves the so-called double-pole approximation (DPA). 
This approximation restricts the complete pole-scheme expansion to the term with the high- 
est degree of resonance. In the case of ly-pair production only the double-pole residues are 
hence considered. The intrinsic error associated with this procedure is aT^/ (nMw) < 0.1%, 
except far off resonance, where the pole-scheme expansion cannot be viewed as an effective 
expansion in powers of Ty /My-, and close to phase-space boundaries, where the DPA cannot be 
trusted to produce the dominant contributions. In the latter situations also the implementa- 
tion dependence of the double-pole residues can lead to enhanced errors. Close to the nominal 
(on-shell) VT-pair threshold, for instance, the intrinsic error is effectively enhanced by a factor 
Mw/{\/s — 2Mw)- In view of this it is wise to apply the DPA only if the energy is several Tw 
above the threshold. 

In the DPA one can identify two types of contributions. One type comprises all diagrams 
that are strictly reducible at both unstable IV-boson lines (see Fig. |l]). These corrections are 
therefore called factorizable and can be attributed unambiguously either to the production 
of the ly-boson pair or to one of the subsequent decays. The second type consists of all 
diagrams in which the production and/or decay subprocesses are not independent and which 
therefore do not seem to have two overall W propagators as factors (see Fig. |]). We refer to 
these effects as non-factorizable corrections.^] In the DPA the non-factorizable corrections arise 
exclusively from the exchange or emission of photons with < 0{Tw)- Hard photons as 

^It should be noted that the exact spHt-up between factorizable and non-factorizable radiative corrections 
requires a precise (gauge-invariant) definition. We will come back to this point in Sect. 
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Figure 2: Examples for virtual (top) and real (bottom) non-factorizable corrections to H^-pair 
production. The black circles denote the lowest-order Green functions for the production of 
the virtual W^-boson pair. 

well as massive-particle exchanges do not lead to double-resonant contributions. The physical 
picture behind all of this is that in the DPA the VT-pair process can be viewed as consisting of 
several subprocesses: the production of the 14^-boson pair, the propagation of the W bosons, 
and the subsequent decay of the unstable W bosons. The production and decay are "hard" 
subprocesses, which occur on a relatively short time interval, 0{1/Mw). They are in general 
distinguishable as they are well separated by a relatively big propagation interval, 0{1/Tw)- 
Consequently, the corresponding amplitudes have certain factorization properties. The same 
holds for the RC to the subprocesses. The only way the various stages can be interconnected is 
via the radiation of soft photons with energy of (9(r^). These photons induce relatively long 
range interactions and thereby allow the various subprocesses to communicate with each other. 

Within the DPA the generic form of the virtual and soft-photonic 0{a) RC to off-shell 
W^-pair production can now be cast in the following gauge-invariant form: 

c^c^DPA = da^pjs^il + (5dpa), (9) 

where (iapp^ is the lowest-order differential cross-section in DPA. The hard-photon effects in 
DPA are added separately, in view of their dependence on the phase-space of the hard photons. 
Their contribution strongly depends on the distribution that is studied, i.e. on the integrations 
that have to be performed. More details will be given in Sect. ^. Since one also knows the exact 
Born cross-sections for off-shell VT-pair production [(ia^^] and for the background diagrams 
contained in the four-fermion process (|) [c^cr^kg] , one can also add those to the above expression. 
The final gauge-invariant result up to 0{a) or, equivalently, 0(Tw/Mw) precision reads 

da = da^pjs^il + 5dpa) + {da^^ - da^pj^) + da^^^^. (10) 

The purpose of this paper is to give a detailed discussion of 6bpa, i-e. the 0{a) corrections to 
dcjpp^. In order to make contact with experimental cross-sections the other terms, da^^^^ — 
dcTppA and da^^^ , are also relevant. The full gauge- invariant Born term, including all diagrams. 
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has been discussed in the hterature |Q. It has also been compared with the non-gauge-invariant 
cross-section da'^^, calculated in the unitary gauge. In many cases cicr^w gives numerically 
a good approximation to the complete Born cross-section. Moreover, in practice it is often 
extracted from the data in the experimental analyses. Therefore it is useful to present a 
numerical comparison between da^^j^^ and daj^-pj^. This will be done in Sect. H. 



3.1 The double-pole approximation: conventions and an example 

As is clear from the above-given discussion of the DPA, a specific prescription has to be given 
for the calculation of the DPA residues. Or, in other words, we have to fix the implementation 
of the mapping of the full off-shell phase space on the kinematically restricted (on-resonance) 
one. We have opted to always extract pure double-pole residues. This means in particular 
that after the integration over decay kinematics and invariant masses has been performed the 
on-shell cross-section should be recovered. 

In the rest of this subsection we will explain our method in more detail by applying it to 
the lowest-order reaction 

e+(gi) e-fe) - W+{p{) W-{p2) - fi{ki)m[) f2{h)K{K), (11) 

involving only those diagrams that contain as factors the Breit-Wigner propagators for the W'^ 
and bosons. Here /i and f[ are the decay products of the boson, and /2 and /2 those 
of the W~ boson. It should be noted that a large part of the RC in DPA to this reaction can 
be treated in a way similar to the lowest-order case, which is therefore a good starting point. 
The amplitude for process (|Tl]) takes the form 



A^^EnWM.M.)^^, (12) 

where any dependence on the helicities of the initial- and final-state fermions has been sup- 
pressed, and 

A = - + iMwTw, = {ki + k'f. (13) 

The quantities A^^\Mi) and A^^\M2) are the off-shell W-decay amplitudes for specific spin- 
polarization states Ai (for the W~^) and A2 (for the W~), with Aj = (— 1,0, +1). The off-shell 
ly-pair production amplitude nAiA2(^i) ^2) depends on the invariant fermion-pair masses Mi 
and on the polarizations Aj of the virtual W bosons. In the limit —>■ Mw the amplitudes 11 
and A^^) go over into the on-shell production and decay amplitudes. 

The choice of the polarization states labelled by Aj is in principle free. The amplitude M. 
is obtained by summing over the polarizations and is therefore independent of such a specific 
choice. As it turns out, it will be convenient to use different choices in different parts of the 
RC calculation (see Apps. |^ and ^). 



In the cross-section the above factorization leads to 



■^1 = L ^[AiA,][A;Ay(Mi,M2) ^ in |2 ■ (14) 



fermion hclicities Ai,A2,A'j^,A2 
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In Eq. (14) the production part is given by a 9 x 9 density matrix 

P[Aa.][A;A^](Mi,M2) = E n,,,,(Mi,M2)ni,,,(Mi,M2). (15) 

helicities 

Similarly the decay part is governed by 3 x 3 density matrices 

^x,K(Mi) = J2 A,, (MO A*,, (MO, (16) 

fermion helicities 

where the summation is performed over the hehcities of the final-state fermions. 



It is clear that Eq. ([15| ) is closely related to the absolute square of the matrix element for 
stable unpolarized VT-pair production. In that case the cross-section contains the trace of the 
above density matrix 

Tr V{Mw,Mw) = J2'P[>^i>^2][XiX2]{Mw,Mw) = E {Hx^x^Mw, Mw)\''. (17) 

Ai,A2 all polarizations 

The decay of an unpolarized on-shell W boson is determined by 

TrV{Mw) = J2^x^xAMw) = E \^xAMw)\'. (18) 

Aj all polarizations 

Note, however, that also the off-diagonal elements ofV{Mw, Mw) and V{Mw) are required for 
determining Eq. ([T^) in the limit Mj — >■ M^. 

As a next step it is useful to describe the kinematics of process (|TTD in a factorized way, 
i.e. using the invariant masses Mi and M2 of the fermion pairs. The differential cross-section 
takes the form 

d^ = Ys^\M? = ^ E l-^r ^rpr ■ dTl, . rfr,-, ■ ^ ■ (19) 

where dT^f indicates the complete four-fermion phase-space factor and s = {qi + q2)'^ the centre- 
of-mass energy squared. The phase-space factors for the production and decay subprocesses, 
dT pj. and dV^^^, read 

■ff-i. = (2^)j*' lit' 

When the factorized form for l-^P is inserted one obtains 

= j; E ^'P[MX,][Kx',]iM,,M2)dTp,xV^^y^iM,)dTl^xV^^^^^^ 

1 dMf 1 dMi ^21) 



2tt |Di|2 2tt 
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As mentioned before, the definition of the DPA residues is not unique. To define them we 
first organize the four-fermion kinematics in a special way. In the laboratory (LAB) frame 
we write the four-fermion phase space in terms of a solid production angle for the W'^ boson 
and sohd decay angles for two of the final-state fermions, one originating from the W'^ boson 
and one from the W~ boson. These angles will be kept fixed at any time during the process 
of defining the DPA residues. For later use we explicitly write down our conventions for the 
momenta, invariants, and phase-space factors. The momenta read 

= £;(!, sin ^,0, cos ^), ^2 = ^(1, - sin ^, 0, - cos^), (22) 
pi = E,{1, 0, 0, ^), p2 = E2{1, 0, 0, -|-), 

-C/l -C/2 

ki = -£3(1, sin 6*3 cos 03, sin 6*3 sin ^3, cos 6*3), /c2 = -^4(1, sin 6*4 cos(j>4, sin 6*4 sin ^4, — cos 6*4), 



with 



E = ly/^, E,,2 = :r^{s + Ml,-Ml,), ^3,4 ^ 



2 2Vs' ' 2 £;i,2 -pcos6'3,4 

p = ^ \^'\^/~s, Ml, M2), A(^, Ml, M2) = [s- (Ml + M^fWs - (Ml - M^f]. (23) 

The momenta of the other final-state particles follow from k[ — pi — ki. The masses of the 
fermions are neglected whenever possible. This, hence, excludes situations in which the fermion 
masses are needed to regularize singularities from the radiation of coUinear photons. The 
Mandelstam invariants are defined as 

s = (gi + ?2)', i = (pi-gi)', u^{j)2-qif^Ml + Ml-s-t. (24) 

Prom all this it should be clear that the invariant masses Mj only appear explicitly in the 
energies and velocities of the W bosons and their decay products. 

In this notation the production phase-space factor reads 

with dTt denoting the solid angle between the boson and the positron. The decay phase- 
space factors are given by 

and a similar expression for dV^^^. Here dfl^ denotes the solid decay angle between the 
boson and /i. 

Our choice of the DPA residues amounts to a two-step procedure for fixing the invariant 
masses Mi and M2, appearing in the four-fermion kinematics and the amplitudes n(Mi,M2), 
A(+)(Mi), and A(-)(M2). The first step is the replacement 

M^^M^- tMwTw, (27) 
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i.e. the residue is taken at the Breit-Wigner poles [see discussion below Eq. (H)]. Note that this 
replacement, of course, does not apply to the Breit-Wigner resonances themselves. The phase- 
space conventions displayed above fix our choice for the implementation. The solid angles are 
kept fixed, whereas the energies and velocities become complex [as can be seen from Eq. (p^)]. 
Note that the so-obtained set of momenta preserves momentum conservation. This protects 
the DPA residues against effectively crossing the four-fermion phase-space boundaries, which 
might lead to a reduced quality of the DPA. 

For practical purposes, however, it is messy to evaluate the 0{a) RC to the amplitudes 
Il{Mi, M2), A^~^\Mi), and A^~\M2) at the complex Breit-Wigner poles. This would require 
the analytic continuation of the one-loop expressions to the second Riemann sheet. As such we 
approximate the DPA residues by using 

Mi = Mw. (28) 

The error introduced by the on-shell approximation (pS]) is of order 0{Tw / ^w)- When this 
error comes on top of the 0{a) RC it can be neglected, since terms of 0{aVw / Mw) are 
neglected anyhow in the DPA approach. By having fixed the solid angles in the mapping, the 
on-shell phase-space points defined with Eq. (|^) remain physical. From this point of view 
it is a sound implementation procedure. A procedure that fixes Mandelstam variables in the 
mapping involves phase-space regions where it may be regarded as being unsound, as has been 
indicated in the example below Eq. (^. In such cases it is preferable to set the cross-section 
to zero rather than evaluating it for unphysical values. 

The thus-obtained amplitudes become (gauge-invariant) on-shell ones, whereas the four- 
fermion phase space is reduced to the phase space of two fermion pairs with invariant masses 
Mt^. Since the DPA forces us to only consider collider energies that are several T\y above the 
on-shell TV-pair threshold, the W-hosow velocities stay well defined in our approximation. The 
only place where the invariant masses Mi still show up is in the Breit-Wigner resonances D,, 
which can be pulled out as overall factors. The integration over the Breit-Wigner resonances 
need not be restricted to the physical region Mj > 0, Mi + M2 < y^. Since the DPA is anyhow 
not valid far off resonance, we can just as well integrate over the full range of the distributions, 
(— cxD, +cxo), which guarantees that the on-shell results are recovered when the decay kinematics 
and invariant masses are integrated out. This means that the following integral will be used: 

— / dMf-\ = ^ . (29) 

27r J ' 2MwTw 

— CO 

So far we have explained how to calculate the canonical multidifferential cross-section 
da/{dfldQ3dflidMfdM2) in the DPA. In case one needs a multidifferential cross-section in 
other variables one should relate that cross-section to the canonical one by means of a Jaco- 
bian. In order to obtain this Jacobian the off-shell kinematical relations (^) and ( ^31) should 
be used. 
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4 A numerical comparison of different Born cross-sections 



In order to have an idea of the differences between the exact Born cross-section da^^^, corre- 
sponding to process (|^), and its DPA limit (iaop^, we now present a brief numerical comparison. 

First we discuss the set of parameters used to produce the plots throughout this paper. To 
facilitate cross-checks with the results presented in the literature, we adopt the LEP2 input- 
parameter scheme advocated in Ref. 0]. In this scheme the Fermi constant G^, the fine- 
structure constant a, and the masses of the light fermions^ and W, Z bosons are the independent 
input parameters. The mass of the top quark, mt, follows from the Standard Model prediction 
for muon decay 

Q = ^ 3_ ^30) 

^ v^M2,(l-M2,/M|)l-Ar- ^ ^ 

The quantity Ar denotes the loop corrections to muon decay. It is zero at tree level, but 
when loop corrections are included it depends on the input parameters as well as on nit, the 
Higgs-boson mass M^, and the strong coupling a^.Q 

In analogy to Ref. we use in our numerical evaluations the following set of (slightly 
outdated) input parameters: 

a = 1/137.0359895, = 1.16639 x 10"^ GeV'^ 

Mz = 91.1884 GeV, Mw = 80.26 GeV, 

me = 0.51099906 MeV, = 105.658389 MeV, = 1.7771 GeV, 

mu = 47 MeV, = 47 MeV, 

m, = 150 MeV, = 1.55 GeV, 

m6 = 4.7GeV, (31) 

and choose 

Mh = 300 GeV, as{Ml) = 0.123. (32) 

The solution of Eq. (|30D, using a calculation of Ar that contains all the known higher-order 
effects, gives the value for the top-quark mass @] 

mt = 165.26 GeV. (33) 



As was already mentioned in Sect. |1|, the use of instead of a in the lowest-order cross- 
sections very often reduces the size of the one-loop non-photonic RC considerably. This so-called 
G^-parametrization consists in the replacement 

da = da\l + (5^-i°°P) ^ _ ^^^^ (1 + - nAri-'°°P) = d(x\l + 5^-i°°p), (34) 

^The masses of the light quarks are adjusted in such a way that the experimentally measured hadronic 
vacuum polarization is reproduced. 

'^The so-obtained top-quark mass will become as- and Af //-dependent. It can be confronted with the direct 
measurements at Fermilab and the indirect ones from LEP in order to obtain indirect information on Mh- 
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where da^ oc a" and according to Eq. (|30|) da^ oc G^. The results presented in this paper are 
all calculated in this parametrization. 

Another important parameter featuring in our calculations is the width of the W boson. 
As explained below Eq. (H), we will use the calculated on-shell width. Since we want the 
Breit-Wigner resonances to be as close to reality as possible, we will always use the 0{a)- 
and C(Q;s)-corrected width Tw, regardless of the fact that we sometimes consider lowest-order 
distributions. Using the above set of input parameters we find in the G^-parametrization 

Ely = 2.08174 GeV. (35) 

For future use we note that the lowest-order VT-boson width in this parametrization reads 
F^ = 2.03540 GeV. It is also relevant to stress that the 0{a) corrections to the leptonic partial 
widths Tw^eui are small and negative (~ —0.3%). The 0{as) corrections to the hadronic partial 
widths are positive, leading to the positive overall correction to the total W^-boson width. 

Having fixed the input, we now compare da^^^ and rfcr^pA for the total cross-section Utot 
(in Fig. ^ and the differential cross-section da/dcos6 (in Fig. where 6 is the polar angle 
between the boson and the positron in the LAB frame [see Eq. (^2]) ]. The latter distribution 
is given for 2E = 184 GeV, whereas cxtot is presented for a range of LEP2 energies. We select 
one particular purely leptonic final state, ^^v^r^Vr- In view of the massless treatment of the 
final-state fermions and the universal lowest-order interaction between the fermions and the 
W bosons, the results for the various final states can be obtained by multiplying the purely 
leptonic result by a factor A*"^^ lKf2/2p- Here Vyj. is the mixing matrix and the 

colour factor. For leptons only Vp^i = 1 is non- vanishing and = 1. 

We consider the following four cases: 



i) The calculation for stable W bosons, multiplied by the branching ratio (r^_^^j^ /F 



WJ 



ii) The DPA calculation, where in Eq. (^Tj) the on-shell approximation is applied to both 
the matrix elements and the four-fermion phase space. The integrations over the 
Breit-Wigner resonances are extended to the full range (— oo, +oo), i.e. Eq. (E^) is used. 



iii) The calculation where the matrix element (|12]) is on-shell, but the four-fermion phase 
space in Eq. (^) is not. The Mf integrations are performed in the physical region. 



iv) The off-shell calculation according to Eq. (PT|), with the M/ integrations performed in the 
physical region. This corresponds to da^^^. 



In cases ii), iii), and iv) the IV-boson propagators contain the width Tw, as given in Eq. (^5]). 
The matrix element in case iv) is not gauge-invariant, it is calculated in the unitary gauge. 

For the total cross-section (see Fig. §) cases i) and ii) differ by a fixed overall factor, as 
was to be expected from the Breit-Wigner integrals. The overall factor is determined by 
{Tw /T^y = 1.04605 and is found to be 1.04609. The agreement with this overall factor is a 
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Figure 3: Comparison of different Born approximations for the total cross-section cxtot as a 
function of the accelerator energy. The four curves correspond to the cases i) - iv) introduced 
in the text. 
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Figure 4: The same curves as in the previous plot, this time for the lowest-order production- 
angle distribution da /dcos 9 at 2E = 184 GeV. 
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check on the numerical integration over the decay angles. For the production-angle distribution 
(see Fig. ^ the same overall factor is observed. 

One of the ingredients of Eq. (p!0| ) is the difference da^/^r — (^Cdpa- By comparing cases ii) 
and iv) in Fig. ^ one observes that this difference varies between 18% and 5%, when going from 
165 GeV to 200 GeV. This is in good agreement with the expected 0{Tw/ AE) precision of 
the DPA limit, where AE is defined as the distance in energy to the VF-pair threshold, AE = 
y^ — 2Mvy. Larger differences arise when the comparison between (Jiyiy and cx^p^ is carried out 
at much higher energies. The reason is that (JiYyy is not SU{2) gauge-invariant and does not 
fall off sufficiently fast. By properly combining cr^^yy and the background contributions a^^^ the 
SU (2) gauge invariance can be restored and cr^pA again turns out to be a good approximation 
to the total (four-fermion) cross-section. 



5 Radiative corrections in the double-pole approxima- 
tion 

In this section we discuss the complete 0{a) RC to process dn]) in the context of the DPA. 
As mentioned in Sect. |^, the most economic way of calculating the RC up to 0{a) precision 
involves the DPA. In this approximation all virtual corrections can be classified into two groups: 
factorizable and non-factorizable corrections. The factorizable corrections are directly linked to 
the electroweak one-loop RC to the on-shell production and on-shell decay of the W bosons. The 
remaining non-factorizable virtual corrections can be viewed as describing interactions between 
different stages of the off-shell process. They will only originate from certain photonic loop 
diagrams, as stated in Sect. ^ The real-photon corrections can also be classified in factorizable 
and non-factorizable corrections, although the various regimes for the photon energy require 
some special attention. 

In the following we will describe the calculation of all these corrections in more detail and 
comment on the accuracy and applicability of the results. At this point we remind the reader 
that throughout the calculations of the RC in the DPA two additional approximations are 
used. First, whenever possible we consider the initial- and final-state fermions to be massless, 
i.e. excluding the cases in which the fermion masses are needed to regularize singularities from 
the radiation of collinear photons. The error of this approximation is at most 0{amT-/Mw) or 
0{a \ Vcb \ rrib/Mw), which is beyond the accuracy of our calculation. Second, we assume that the 
accelerator energy is sufficiently far (read: several Tw) above the threshold for on-shell VT-pair 
production. Close to threshold the DPA cannot be trusted to produce the dominant contri- 
butions and therefore our approach breaks down. The accuracy of this "far-from-threshold" 
approximation is 0{arw/AE), where AE is the distance in energy to the TV-pair threshold, 
AE= y/^- 2Mw. 
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5.1 Virtual corrections 



As a first step we discuss how to separate the virtual corrections into a sum of factorizable 
and non-factorizable virtual corrections. The diagrammatic split-up according to reducible and 
irreducible VF-boson lines is an illustrative way of understanding the different nature of the 
two classes of corrections, but since the doublc-rcsonant diagrams are not gauge-invariant by 
themselves the precise split-up needs to be defined properly. 

Wc can make use of the fact that there are effectively two scales in the problem: Mw and 
Tw Let us now consider virtual corrections coming from photons with different energies: 

• soft photons, <^ F^y, 

• semi-soft photons, Ej — 0{rw), 

• hard photons, Tw <^ E^ = 0{Mw)- 

Only soft and semi-soft photons contribute to both factorizable and non-factorizable corrections. 
The latter being defined to describe interactions between different stages of the off-shell process. 
The reason for this is that only these photons can induce relatively long-range interactions 
and thereby allow the various subprocesses, which are separated by a propagation interval 
of 0{1/Tw): to communicate with each other. Virtual corrections involving the exchange of 
hard photons or massive particles contribute exclusively to the factorizable corrections. In 
view of the short range of the interactions induced by these particles, their contribution to the 
non-factorizable corrections are suppressed by at least 0{Tw / M^). 

As hard photons contribute to the factorizable corrections only, we only need to define a 
split-up for soft and semi-soft photons. It is impossible to do this in a consistent gauge-invariant 
way on the basis of diagrams. As we will see, it might happen that only part of some particular 
diagram should be attributed to the non-factorizable corrections, the rest being of factorizable 
nature. 

The matrix element for soft and semi-soft photons can be written as a product of the 
lowest-order matrix element in DPA (A^dpa) and conserved (semi) soft-photon currents. These 
currents can be decomposed into production and decay currents with the help of a partial- 
fraction decomposition for virtual-photon emission from a VF-boson line: 

[p2 - M2] [p2 + 2pk - M2 + io\ ~ 2pk + io [y^^AP ~ p^ + 2pk -M'^ + io 



2pk + io \D D + 2pk + io 

Here = — iMw^w, k is the loop momentum of the exchanged photon, p + k is the 
momentum of the W boson inside the integral, and p is the momentum of the W boson outside 
the integral. The infinitesimal imaginary part +io is needed to ensure a proper incorporation 
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Figure 5: Examples of one- loop diagrams that contribute to both factorizable and non- 
factorizable corrections. 



of causality. In this way one obtains a sum of two resonant VT-boson propagators multiplied by 
an ordinary on-shell eikonal factor. This decomposition allows a gauge-invariant split-up of the 
matrix element in terms of one contribution where the photon is effectively emitted from the 
production part, and two others where the photon is effectively emitted from one of the two 
decay parts. The squares of the three contributions can be identified as factorizable corrections, 
whereas the interference terms are of non-factorizable nature. 

In order to illustrate our method, we explicitly apply it to two special one-loop contributions 
(see Fig. The first one [diagram (a)] is the so-called Coulomb graph, involving photon 
exchange between the two W bosons. The corresponding semi-soft matrix element reads 



DPA 



DPA 



(27r)4[A;2 + io] [Di + 2kpi + io] [D2 - 2kp2 + io] 

d'^k 4:pip2 \ Di 



(27r)4[A;2 + io] [2kpi + io] [-2kp2 + io] { Di + 2kpi + io 



D2 D1D2 



D2 - 2kp2 + 10 [Di + 2kpi + io] [D2 - 2kp2 + io]}' 



(37) 



The first term in the last expression gives rise to a factorizable (on-shell) contribution to the 
production stage, whereas the other three terms are counted as non-factorizable contributions. 
These three terms are classified as (prod x dec"*"), (prod x dec"), and (dec"*" x dec"), respectively. 
In other words, the Coulomb graph contributes both to the usual Coulomb effect in on-shell 
W^-pair production and to a non-factorizable part. For the photonic interaction between the 
boson and its decay product [diagram (b)] we obtain 



Mb = -ie^M^opj^ J 



d^k 4pifci 



(27r)4[A;2 + io] [2kki + io][Di + 2kpx + io] 



-/ ^^^j^2'KY[k^ + io][2kpi+io][2kki+io]\ Di + 2kpi + io]' ^ ' 

Again the first term is a factorizable contribution, belonging to the decay stage, whereas 
the second term is a non-factorizable contribution of the type (prod x dec"*"). As is clear 
from these examples, the non-factorizable contributions always involve the Breit-Wigner ratios 
Di/{Di ± 2kpi), which effectively remove the overall W propagator 1/Di. The more energetic 



18 



the exchanged photon is, the more suppressed such a ratio will be in the vicinity of the Mp reso- 
nance. In fact, for > A the non-factorizable contributions are suppressed by 0{Mwrw/[EA]) 
(see App. 0). 

Every one-loop diagram with a semi-soft photon can be treated in this way. By collecting 
all terms that contain the ratios Di/[Di ± 2kpi] the formula for non-factorizable corrections 
is obtained. As one can see explicitly below, this expression is gauge-invariant. Since the 
expression contains those diagrams where irreducible VF-boson lines are present, it can be 
viewed as a gauge-invariant extension of the set of IV-irreducible diagrams. 



The so-defined non-factorizable corrections read 

d^k 



^^"^ ' (27r)4[fc2 + to] 
The currents are given by 



(39) 







Pi 



+ 



1^ 

P2 



?2 



kqi + io kq2 + io 



kpi + io —kp2 + io 



9i 



(l2 



kqi + io —kq2 + io 



(40) 



for photon emission from the production stage of the process, and 



kpi + io 
-kp2 + io 



+ Qh 



kki + io 



k[ 



kk[ + io 



+ Qh 



-kk2 + io 



+ 2kpi ' 

D2 



—kk'n + io 



D2 - 2kp2 



(41) 



for photon emission from the decay stages of the process. Here Qj stands for the charge of 
fermion / in units of e. After having defined the gauge-invariant currents, the +io terms can 
be dropped from Di + 2kpi and D2 — 2kp2, since ImDi > 0. Note the difference in the sign 
of the io parts appearing in the currents JTq and J'^. These signs actually determine which 
interference terms give rise to a non-vanishing non-factorizable contribution after virtual and 
real-photon corrections have been added. As can be seen from Eq. (|38|), in the upper hemisphere 
of the complex ko plane there is only one pole: the so-called 'photon pole', originating from the 
photon propagator l/lk"^ + io]. When virtual and real-photon corrections are combined, such 



a non-factorizable contribution will vanish |21|. For the Coulomb graph, Eq. (37), this is not 



the also poles from the other propagators are present in both hemispheres. As a result 

of such considerations only a very limited subset of 'final-state' interferences survives [^, p!0| : 
the virtual corrections corresponding to Figs. ^ and ^(a) as well as the associated real-photon 
corrections. 

The virtual factorizable corrections consist of all hard contributions and the above-indicated 
part of the semi-soft ones. The so-defined factorizable corrections have the nice feature that 
they can be expressed in terms of corrections to on-shell subprocesses, i.e. the production of 
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two on-shell W bosons and their subsequent on-shell decays. The corresponding matrix element 
can be expressed in the same way as described in Sect. ^ 

XSS^En.UM.M.)^^. (42) 

Here two of the amphtudes are taken at lowest order, whereas the remaining one contains all 
possible one-loop contributions, including the W wave-function factors that appear in Eq. (|]). 
In this way the well-known on-shell RC to the production and decay of pairs of W bosons 



[[15|] appear as basic building blocks of the factorizable corrections. Q In the semi-soft limit the 
virtual factorizable corrections to the production stage, contained in 11, will cancel against the 
corresponding real-photon corrections. Non-vanishing contributions from 11 occur as soon as 
the k"^ terms in the propagators cannot be neglected anymore. An example of this is the fac- 
torizable correction from the Coulomb graph, given in Eq. (pT]). For the on-shell (factorizable) 
part of the Coulomb effect photons with momenta ko = 0{AE) and \k\ = 0{\/Mw^E ) are 
important |Tl[, i.e. fc^ cannot be neglected in the propagators of the unstable particles. Since 



we stay well away from the ly-pair threshold {AE T\y), this situation occurs outside the 
realm of the semi-soft photons. This fits nicely into the picture of the production stage be- 
ing a "hard" subprocess, governed by relatively short time scales as compared with the much 
longer time scales required for the non-factorizable corrections, which interconnect the different 
subprocesses. 



5.2 Real-photon radiation 

In this subsection we discuss the aspects of real-photon radiation in the DPA. To this end we 
consider the process 

e+(gi)e-(g2) - W+{p^)W-{p2) b{k)] ^ Uk,)fl{k',) h{k2)m,)^{k), (43) 

where in the intermediate state there may or may not be a photon. We will show how to 
extract the gauge-invariant double-pole residues in different situations. The exact cross-section 
for process (H3) can be written in the following form 



2Re{MoM\+MoM*_+M+M*_] + \MoY + \M^Y + \M 



(44) 

where dT^f^ indicates the complete five-particle phase-space factor, and the matrix elements 
Aio and A4± correspond to the diagrams where the photon is attached to the production or 
decay stage of the three W^-pair diagrams, respectively. This split-up can be achieved with the 
help of the real-photon version of the partial- fraction decomposition (|5B|). Each contribution 



^Note that the complete density matrix is required in this case, in contrast to the pure on-shell calculation 
which involves the trace of the density matrix only. 
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to the cross-section can be written in terms of polarization density matrices, which originate 
from the amphtudes 

^„.n,(M.A«^^. ,45) 
^,.n,M,,.M,^»^. ,46) 



A(+)(Mi) 


A(-)(M2) 








1 A(-)(M2) 




D2 


AW(Mi) 






D2^ 



= n(Mi, M2,) „ " , (47) 

where all polarization indices for the W bosons and the photon have been suppressed, and 

A7 = A + 2kki + 2kk'„ Ml = Mf + 2kki + 2kk[, = {h + k[)\ (48) 

The matrix elements 11^ and A^^^ describe the production and decay of the W bosons ac- 
companied by the radiation of a photon. The matrix elements without subscript 7 have been 
introduced in Eq. (|l^). 

In the calculation of the Born matrix element and virtual corrections only two poles could 
be identified in the amplitudes, originating from the Breit-Wigner propagators The 
pole-scheme expansion was performed around these two poles. In contrast, the bremsstrahlung 
matrix element has four in general different poles, originating from the four Breit-Wigner 
propagators 1/A and 1/Di^. As mentioned above, the matrix element can be rewritten as a 
sum of three matrix elements {M.q, M.j^, M.-), each of which only contain two Breit-Wigner 
propagators. For these three individual matrix elements the pole-scheme expansion is fixed, 
as before, to an expansion around the corresponding two poles. However, when calculating 
cross-sections [see Eq. (0)] the mapping of the five-particle phase space introduces a new type 
of ambiguity. The interference terms in Eq. (^^ involve two different double-pole expansions 
simultaneously. As such there is no natural choice for the phase-space mapping in those cases. 
The resulting ambiguity (implementation dependence) might have important repercussions on 
the quality of the DPA calculation and therefore deserves some special attention. In this context 
the three earlier-defined regimes for the photon energy play a role: 



• for hard photons [E^ ^ '^w\ the Breit-Wigner poles of the W-boson resonances before 
and after photon radiation are well separated in phase space (see Mf^ and defined 
above). As a result, the interference terms in Eq. (|^) can be neglected. This leads 
to three distinct regions of on-shell contributions, where the photon can be assigned 
unambiguously to the W-pair-production subprocess or to one of the two decays. This 
assignment is determined by the pair of invariant masses (out of Mf and Mf^) that is 
m the region. Therefore, the double-pole residue can be expressed as the sum of 
the three on-shell contributions without increasing the intrinsic error of the DPA. Note 
that in the same way it is also possible to experimentally assign the photon to one of the 
subprocesses, since misassignment errors are suppressed, assuming for convenience that 
all final-state momenta can ideally be measured. 
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• for semi-soft photons [E^ = OlT^/)] the Breit-Wigner poles are relatively close together 
in phase space, resulting in a substantial overlap of the line shapes. The assignment of the 
photon is now subject to larger errors. Moreover, since the interference terms in Eq. 
cannot be neglected, the issue of the phase-space mapping has to be addressed. In the 
following we give a proper prescription for calculating the DPA residues and discuss their 
quality. 

• for soft photons [E^ <^ Tw] the Breit-Wigner poles are on top of each other, resulting in 
a pole-scheme expansion that is identical to the one without the photon. 



5.2.1 Hard photons 



Let us first consider the hard-photon regime in more detail. Due to the fact that the poles are 
well separated in the hard-photon regime, it is clear that the interference terms are suppressed 
and can be neglected: 



da = — 
2s 



dT 



4/7- 



(49) 



Note that each of the three terms has two poles, originating from two resonant propagators. 
However, these poles are different for different terms. The phase-space factor can be rewritten 
in three equivalent ways. The first is 



dr 



4/7 



dvl 



dV^j. ■ dV^^^ ■ dV^^^ 



dMl dMl 



27r 



27r 



with 



(27r 



1 rfpi dp2 

o[qi + q2-pi-p2- k) 



dk 



2pw 2p2o (27r)32fco 



The two others are 



dr 



4/7 



dV 



+ 



dTpj. ■ dT^^^ ■ dr^^^ 



dMf^ 

271 



dMl 

271 ' 



with 



dec 



1 



6{pi — ki — k[ — k) 



dki dk[ dk 



(50) 



(51) 



(52) 



(53) 



(27r)2 ^ i ^ 2kio 2k[Q (27r)32fco 

and a similar expression for dFl. The phase-space factors dVpj- and dV^^^ are just the lowest- 
order ones. The cross-section can then be written in the following equivalent form 



da = — 
2s 



lA^ol' dVl + \M+\^ dTl + \M-\'' dT2 



(54) 



In order to extract gauge-invariant quantities, the DPA limit should be taken. This amounts to 
taking the limit using a prescription that resembles the one presented in Sect. |^. 

Note however that pi^2 can be different according to the 5-functions in the decay parts of 
the different phase-space factors. To be specific, the production term |A^oP has poles at 
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p2 = Mf = M^, \M+\^ has poles at pj = Mf^ = and pi = M| = M^, and |A^_p has 
poles at pI = Ml = and p| = = M^. Again we fix solid angles in the mapping, 
including the solid angle of the photon. Since the energy range of the photon in the on-shell 
kinematics of the W bosons is different from the off-shell case, it may happen that the energy of 
the photon in an off-shell four-fermion - one-photon event with certain solid angles lies outside 
the on-shell phase space with the same solid angles. A possible procedure is to assign a zero 
cross-section to those events. Since the events are anyhow rare, being close to the edge of 
the off-shell phase space, this procedure constitutes an acceptable and simple solution. An 
alternative way to avoid unphysical on-shell phase-space points would be to write the photon 
energy as a fraction of the maximum attainable photon energy for given invariant masses yp? 
of the resonant W bosons and given solid angles: 

E, = , angles ). (55) 

In this way the photon energy is projected on the interval [0, 1]. The maximum photon energy 
^max (jgpg]2fis qu the spccific term in Eq. (|5^). Subsequently the fraction and the afore- 
mentioned solid angles are kept fixed during the mapping from off- to on-shell events. Then 
for the on-shell case is found from Eq. (|55|), where pf are replaced by M^. 



It should be stressed that in the on-shell phase space there is no ambiguity concerning 
the treatment of the photon. One obtains in the DPA limit three gauge-invariant on-shell 
contributions to Eq. (0). The calculation of the corresponding matrix elements is presented 
in App. p. 



5.2.2 Semi-soft and soft photons 



We complete our survey of the different photon-energy regimes by considering semi-soft and 
soft photons. The split-up of factorizable and non-factorizable real-photon corrections proceeds 
in the same way as described in the previous subsection for virtual corrections. The result reads 
in semi-soft approximation 



1 



da = —\M^\'^dT4{^ ^ -da^ 



dk 



2s'- ""°P^(27r)32fco 
The gauge-invariant currents Xq and X± are given by 



2Re(J^j:,^ 



-^0 



Pi 



P2 



kpi kp2 



kqi kq2 



— e 



+ e 



kpi^^'kk, ^^^kk[ 



P2 

kp2 



k^ 



Di + 2kpi 
D2 + 2kp2 ' 



(56) 



(57) 



The first three interference terms in Eq. (^B|) correspond to the real non-factorizable corrections. 
The last three squared terms in Eq. ( p6D belong to the factorizable real-photon corrections. 
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They constitute the semi-soft hmit of Eq. (|54D . Note also that the currents are the same for 
both possible expressions for pj, i.e. it does not matter whether pi = ki + k[ or pi 

As mentioned before, the DPA residues have to be defined properly in the presence of 
overlapping Breit-Wigner resonances in the semi-soft regime. The above equation (^) specifies 
such a procedure: the Born cross-section for a four-particle phase space is factored out and 
does not depend on the photon momentum. The factor between the square brackets is the 
usual soft-photon factor, except that the rapid variation of the Breit-Wigner resonances has 
been kept, leading to the ratios Di/Di^ which take into account the shift in the Breit-Wigner 
distributions due to the photon. In the vicinity of the Mf resonance these ratios are negligible 
for hard photons, unity for soft photons, and of 0{1) for semi-soft photons. 

This prescription is by no means unique. In principle one could have chosen any of the two 
Breit-Wigner resonances, 1/Di or 1/Di^, to define the phase-space mapping for the interfer- 
ence terms. Or in other words, both on-shell phase-space limits, Mf = or M^^ = My^, 
constitute equally plausible DPA mappings. The differences between both phase-space map- 
pings are of 0{k). Since the interference terms are only non-negligible in the semi-soft regime, 
it is conceivable that the implementation dependence associated with these different on-shell 
limits remains of 0{rw/Mw)- Let us, for instance, consider the Mf distribution in the vicin- 
ity of the resonance. Any 0{k) shift in the factor multiplying 1/Di would result in additional 
terms of 0{k/Mw), i-e. at worst 0{rw /Mw) in the semi-soft limit. A similar shift in the factor 
belonging to 1/ Di^ results in additional contributions to the DPA residues that are suppressed 
by 0{rw/E) for all values of the photon energy, as is explained in App. 0. Therefore one can 
conclude that overlapping Breit-Wigner resonances do not necessarily imply a reduced quality 
of the DPA. 

Based on this observation, we have the freedom to choose a suitable procedure for semi-soft 
photons. The fact whether hard-photon contributions are suppressed or not will serve as our 
guideline for fixing the choice. Whenever the hard-photon effects yield vanishing contributions 
to certain observables, we will use the semi-soft currents of Eqs. (|56|) and (ISTf) . As a consequence, 
the non-factorizable corrections are always calculated with the help of the semi-soft interference 
terms. In order to specify our approach, we indicate in the following how the above-defined 
matrix elements and currents are used in the various distributions of experimental interest. 

5.2.3 Real-photon contributions to various distributions 

As in all RC, the role of the bremsstrahlung process is twofold. In the first place the off-shell 
ly-pair process (^) has at lowest-order level all kinds of distributions to which one would like to 
calculate 0{a) RC. These RC consist of virtual and real-photon contributions. An integration 
over all allowed photon energies should be carried out, i.e. the radiated photon should be treated 
inclusively. 

The second application of the bremsstrahlung process involves the evaluation of exclusive 
photon distributions. Since the photon has to be detected it should be sufficiently hard. How 
hard depends on the experimental resolution. An example of such an exclusive photon distri- 
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^ 
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^ 

interference ([56|) 
soft interference (|5BD 



Table 1: Formulae for the bremsstrahlung matrix elements for the distribution da/dMfdM,^ 



bution is the photon-energy spectrum da/dE^. This distribution receives contributions from 
the three hard-photon terms in Eq. as well as from the semi-soft interference terms of 
Eq. (|56[) . The latter terms of course only contribute for photon energies that are not too hard. 



In the next section results will be given for various distributions. For exclusive photons we 
present the photon-energy spectrum da/dE^. For the inclusive photon distributions we discuss 
in the following how the bremsstrahlung part is treated. Depending on the distribution, different 
approximations can be made. 

As indicated above, the calculation of RC to distributions of the off-shell TV-pair process 
(^) involves an integration over the photon phase space that is left available when fixing other 
kinematic variables. This means that the photon phase space for hard, semi-soft, and soft 
photons will in general be integrated over. The soft-photon part is standard and should be 
combined with the virtual corrections. How the non-soft part should be treated depends on the 
distribution one likes to study. In general one considers distributions of the form da^pA/d^ , 
where dX stands for dMldM"^ and a product of solid angles for W-hoson production and decay. 
For the following discussion it does not matter whether these solid angles are kept or integrated 
out. However, it does matter whether one wants a double distribution in Mf and M|, a single 
distribution in M^, or no distribution in Mf at all. In other words, the treatment differs 
depending on whether one integrates over Mf or not. 

Let us first discuss the procedure for the real-photon corrections to the double Breit-Wigner 
distribution da/dMfdM2 in the vicinity of the peak, i.e. one is not interested in the tails of this 
peaked distribution. In Table |I| we specify which expressions are used in the different regimes 
for the photon energy. For this specification we use the positions of the Breit-Wigner poles as 
listed below Eq. 



The entries of Table |1] can be explained as follows. The real-photon corrections to the 
distribution da/dM^dM2 come in the first place from the production part, where the hard- 
photon matrix element squared |A^oP contains the overall Breit-Wigner line shapes in and 
M|. This term should be taken from Eq. (5|) and should be integrated over the photon phase 



space. In the soft-photon limit the usual soft-photon factorization in terms of |2^| is obtained. 
This explains the second column in Table. |I|, where the usual RC to stable W^-pair production 
appear, except that one now calculates the full density matrices rather than merely the trace. 
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Table 2: Formulae for the bremsstrahlung matrix elements for the distribution da/dMf. 

The second type of contribution involves the VT-boson decays with additional photon. When 
the photon is hard these contributions tend to zero as F^/Zcq (see App. |C|), since they do not 
contain a double Breit-Wigner distribution in both and M| according to the list below 
Eq. (0). For (semi)soft photons, however, one gets overlapping Breit-Wigner resonances and 
the |X^| terms of Eq. (|56D can be used. These (semi) soft-photon factors have to be integrated 
while keeping Di fixed, i.e. we keep the variables Mf fixed. The corresponding integrals can 
be found in App. y. The integration over the semi-soft photon momenta leads to contribu- 
tions from values that are higher then M?, resulting in a distortion of the original Mf 
Breit-Wigner distribution. This final-state-radiation effect has been discussed recently in the 
literature and turns out to be quite sizeable The reason why the distortion is so large 
lies in the fact that the final-state collinear singularities [oc ^ Q'j log^nij / M^y)] do not vanish, 
even not for fully inclusive photons. After all, a fixed value of makes it impossible to sum 
over all degenerate final states by a mere integration over the photon momentum.0 On top of 
that the IR poles oc induce a further enhancement of the collinear logarithms by a factor 
log(A/M2.). 

So far we have included all factorizable corrections. To this one should add the non- 
factorizable corrections. Again effectively only the (semi) soft region is relevant. For more 
information on these non-factorizable corrections we refer to the detailed discussion in the 



literature R 10, 22 . 



As the next relevant distribution we discuss the real-photon corrections to the single Breit- 
Wigner distribution da/dMf, obtained from the previous case by integrating over M|. To 
this end we make a similar table (Table 0) and discuss the necessary changes. Based on the 
discussion below Eq. (0), it is clear that a Breit-Wigner distribution in Mf is present explicitly 
in |A^oP Ei'iid |A^_p, both for hard and semi-soft photons. These two terms should be taken 
from Eq. (|5^) and should be integrated over the photon phase space and M|. The |A^+P term 
does not have a Breit-Wigner distribution in as long as the photon is hard. Therefore it 
effectively only contributes in the (semi)soft regime, like in the previous case. Therefore the 
|X^| term in Eq. (|56D is used and the photon integration is performed while keeping Mf fixed. 
The non-factorizable corrections contribute in the same way as described for the previous case. 

We conclude by considering pure angular distributions. In this case the picture is simple. 



^The usual cancellation of final-state collinear singularities is achieved only if also is integrated over. 
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Figure 6: The photon-energy distribution da/dE^ for the ii^i^^T Ut'J final state at 2E = 
184 GeV. In addition the separate production and decay contributions are given. 



All contributions should be taken from Eq. (0) and should be integrated over the photon phase 
space and the invariant masses Mf. The non-fact orizable corrections vanish in this situation, 
which is typical for these interference effects pT . 



6 Numerical results 



All the corrections discussed in the previous sections were combined into a Fortran program. 
In this section we present some numerical results. In particular we study the implications of 
the RC on various distributions. For the numerical evaluations we use the G^-parametrization 
(see Sect. |) and the input parameters listed in Eqs. (^T])-(^). 



6.1 One-dimensional distributions 

We start off by considering various one-dimensional distributions. Where applicable both the 
lowest-order distribution {da^ppj and the corrected one ((icTDPA) will be presented, followed by 
a separate plot for the relative correction factor 5dpa [see Eq. @]. 
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Figure 7: Decay contributions to the photon spectrum from different leptonic W decays. 
6.1.1 The photon-energy spectrum 

Since real-photon bremsstrahlung contributes to the RC to various distributions, it is useful to 
evaluate the photon-energy spectrum da/dE^ separately as well. In DPA it gets contributions 
from the three terms in Eq. (|5^) and the semi-soft interference terms of Eq. ( ^6]) . The photon 
spectrum originating from the production stage is the same for all final states, but the spectra 
from the decay stages depend on the specific final state. In Fig. ^ the DPA photon-energy 
distribution da/dE^ is shown for the /i+i/^r~z/T-7 final state at 2E = 184 GeV, together with 
the production and decay parts of the spectrum. In Fig. |^ we separately list the three possible 
leptonic radiative-decay contributions, originating from W iuel = e, fi, r). The substantial 
differences are caused by the explicit fermion-mass dependence for coUinear photon radiation 
(see App. p^ ). 



6.1.2 The total cross-section as a function of energy 

In Fig. 1^ we compare the total cross-section with and without RC for the ix'^u^t^Ut final state. 
The corresponding relative correction factor 5 is given in Fig. ^. As a check of our calculation 
we have verified that the corrected cross-section coincides within the integration errors with 
the corrected cross-section for stable W bosons multiplied by the corrected branching ratio 

The observed RC are large and negative, especially close to the IV-pair threshold. This 
is mainly caused by real-photon ISR, which effectively lowers the available ly-pair energy. 
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Figure 8: The energy dependence of the total cross-section atot for the ^'^I'^r Uj. final state. 
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Figure 10: The production-angle distribution da/dcos6 for the fi^u^r final state at 2E = 
184 GeV. 

combined with the fact that near the H^-pair threshold the cross-section is rapidly decreasing 
with decreasing energy. 



6.1.3 Production-angle distribution 

In Fig. ^ we plot the production- angle distribution da/dcos6 for the yU^z/^r"//,- final state at 



2E = 184 GeV. The relative correction factor displayed in Fig. ^ differs substantially from 
the —12.8% normalization effect that was observed for the total cross-section. In the forward 
direction it is slightly more negative, whereas in the backward direction it is substantially less 
negative. A proper understanding of this distortion is quite important, since non-standard triple 
gauge-boson couplings might lead to exactly this type of signature. The origin of the distortion 
can be traced back to hard-photon ISR. Such hard-photon emissions boost the centre-of-mass 
(CM) system of the W bosons, causing a migration of events to take place from regions with 
large cross-sections in the CM system (forward direction) to regions with small cross-sections 
in the LAB system (backward direction). The more peaked the distribution is, the stronger the 
boost effects will be. 



6.1.4 Invariant-mass distribution 



In Fig. |12| we compare the Breit-Wigner line shape da/dMf with and without RC for the 
li^u^T'Vr final state at 2E = 184 GeV. Since the corrected line shape receives completely dif- 
ferent contributions from the production and decay stages, these parts are displayed separately. 
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Figure 11: The relative correction factor corresponding to Fig. |T0 
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Figure 12: The invariant-mass distribution da/dMf for the yU+z/^r Uj. final state at 2E 
184 GeV. 
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Figure 13: The relative correction factors S corresponding to Fig. 12 



The corresponding relative correction factors are presented in Fig. 13. 



The correction to the production stage leads to a constant reduction of the Breit-Wigner 
line shape. The corrections to the decay stages comprise the factorizable decay corrections 
(columns 3 and 4 in Table 0) and the non-factorizable corrections (column 5 in Table 0). The 
latter are very small ||, |10|- We see that also the decay corrections reduce the Breit-Wigner 
line shape. The amount of reduction depends on the particular final state, as can be seen 
from Figs. [1^ and |T5|, where we consider different leptonic final states. The differences are 
caused by the explicit fermion-mass dependence for coUinear photon radiation (see Apps. |B.4| 
and y). For other final states, involving quarks, the reduction is determined by the masses and 
charges of the decay products of the boson (see App. |CD. Besides the reduction, the decay 
corrections also distort the resonance shape, as is clear from Fig. O. This final-state-dependent 



FSR distortion effect has recently been discussed in the literature [0]. It can be quantified in 
terms of the shift in the peak position of the Breit-Wigner line shape, which we find to be —20, 
—39, and — 77MeV for l^+-boson decays into r+z/^, fi'^i'^, and e"'"z/e, respectively. This is in 
agreement with the (leading) shifts predicted by the H^-boson version of Eq. (19) in Ref. 0, 
taking into account the fact that the correction to the production stage reduces the line shape 
by 12.0%. The observed shifts differ by —5 MeV to —10 MeV from these predictions as a result 
of the non-leading terms that are present in our full DPA calculation. 

It should be stressed that theoretically the large distortion is a genuine effect. It would 
also be as relevant in practice if the invariant mass of the fermions could be measured. For 
various reasons this is problematic. One reason is the almost unavoidable inclusion of a collinear 
photon in the measured invariant mass. Such an inclusion would effectively decrease the leading 
logarithm [log(M^/mj)] that dominates the effect. What remains of the distortion in practice 
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Figure 14: Distortion of the invariant-mass distribution da/dMf at 2E = 184 GeV due to the 
decay corrections. Three leptonic final states are considered: r^UrT^i/r (Tr-decay), /i+z/^r"//,- 
(yLir-decay) , and e"'"z/er~i/^ (er-decay). 
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Figure 15: The relative correction factors 6 corresponding to Fig. |14 
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Figure 16: The decay-angle distribution da /d cos 63 for the ^'^i^^t final state at IE = 
184 GeV. Here 63 = l{iJ,+, W+) in the LAB frame. 



should be studied with an event generator, which simulates the actual measurement. 



6.1.5 Decay-angle distribution 



In Figs. and the results are shown for the decay- angle distribution da /d cos 63 for the 
fi'^Uf^T^UT- final state at 2E = 184 GeV. Here 63 is the angle between the /i"*" and the W'^ boson 
in the LAB frame [see Eq. (P^J. The correction is negative and becomes more negative for 
forward angles. The shape of the relative correction factor is again the result of hard-photon 
boost effects. 



6.2 Double invariant-mass distribution 



Finally we consider the two-dimensional distribution da/dM^dM^ at 2E = 184 GeV, evaluated 
using the contributions specified in Table 0. Instead of plotting the absolute distributions, only 
the relative correction factors S are presented. We do this for two specific leptonic final states. 
In Fig. |18| the e~^i>ee~Ue final state is considered. This final state has the largest correction. 
Keeping one fixed, the correction to the other distribution is qualitatively the same as the 



correction to the one-dimensional distribution in Fig. For the e~^i'e^~i/T final state (Fig. [ToD 
the correction is clearly not symmetric in the Mf. This was to be expected, since the decay 
corrections for e~*"i^e and T~Ur differ appreciably. 

In Tables ^ and ^ we present some explicit values for the relative correction factor, split up 
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Figure 18: Correction to the double invariant-mass distribution da/dMfdM2 for the e"'"z/ee Uf, 
final state at 2E = 184 GeV. 
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Figure 19: Correction to the double invariant-mass distribution da/dMldM"^ for the e"'"i/eT Vt 
final state at 2E = 184 GeV. 
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Table 3: Relative correction factors [in %] for the double invariant- mass distribution 
da/dMldMl at 2E = 184 GeV. Left: the corrections fr om the W^+-boson decay stage 6'^^^{Mi) 
for different leptonic decay channels and three near- resonant invariant masses. Right: the 
non-factorizable corrections 5ni{Mi, M2). 
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Table 4: Relative correction factors [in %] for the double invariant-mass distribution 
da/dMldM^ at 2E = 184 GeV. Left: the e^Uee'Ue final state. Right: the e+i/gT'^^r final 
state. 

into the separate contributions according to 

The correction from the production stage is constant as a function of the invariant masses, 
5pr = —12.0%. The non-factorizable contribution Snf{Mi, M2) is given in Table ^for three near- 
resonant values for the invariant masses Mj. We indicate these three values by Aj = (Mj — 
Mw)/T]y = —1/2, 0, 1/2. The non-factorizable corrections do not depend on the particular 
leptonic final state. The corrections S^^^{Mi) from the decay stages do depend on the particular 
leptonic final state as explained before. In Table ^ we present these corrections for the three 
leptonic decay modes and the three near-resonant values for the invariant mass. 



7 Discussion and conclusions 

In this paper 0{a) radiative corrections (RC) to four-fermion production in e~^e~ collisions 
have been studied. The energy region chosen is that of LEP2, where the four-fermion final 
state is predominantly formed through intermediate W^-pair production. 

Since a complete 0{a) RC calculation for a two-particle to four-particle process is beyond 
present possibilities, a sensible approximation scheme has to be used. The smallness of Tw/Mw 
offers the possibility to use the double-pole approximation (DPA). In practice it means that we 
calculate 0{a) and 0(Tw/Mw) corrections but neglect 0{aTw/Mw) corrections. Although 
this approximation has been advocated in the literature for quite some time, so far no complete 
0{a) results have been given. As far as we know, the study that has come closest to achieving 
a complete DPA calculation for an actual process involves the factorizable QCD corrections to 
the process e~^e~ — > ti [^. Many of the issues discussed in the present paper have, however, 
not been addressed in Ref. 



We have applied the method to VT-pair-mediated four-fermion production for a number of 
reasons. There is the methodological aspect of dealing with unstable particles in DPA, involving 
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issues like gauge invariance, interactions between different stages of the reaction, RC to density 
matrices, and the phase-space mappings. All of these issues play a role in the ly-pair-mediated 
four-fermion production process. There is also the practical aspect of assessing how large 0{a) 
corrections can be for certain distributions. This is of importance for W^-pair studies at LEP2. 

On the methodological side, we have succeeded in finding a consistent prescription for 
applying the idea of the DPA. The kinematics for calculating the poles of the matrix elements 
is necessarily on-shell kinematics. Also the phase-space factor in the cross-section is treated in 
on-shell kinematics. The off-shell invariant masses occur only in the Breit-Wigner factors. All 
of this is well defined both for the radiative and non-radiative parts of the cross-section and 
therefore our calculation itself is unambiguous. 

An unavoidable problem is the relation between off-shell four-fermion events and the on-shell 
calculated events. This question arises when one likes to connect experimental cross-sections 
to the calculated DPA cross-sections. Also here recipes are chosen both for the radiative and 
non-radiative phase-space points. For the latter the mapping is natural if one chooses the 
invariant masses Mi^2 and angles as variables. All off-shell points can be mapped onto on-shell 
points. For the radiative events in an off-shell phase space, the photon variables have to be 
added. A natural choice is the photon energy and angles. If the mapping is chosen such that 
the photon variables remain the same, one sometimes maps outside the on-shell phase space. 
Different remedies for this problem are possible. One can choose a procedure that assigns zero 
cross-sections to those points. On the other hand, the photon energy in the off-shell four- 
fermion - one-photon phase space may be rescaled in order to obtain a physical point in the 
on-shell phase space. In general, there will be a dependence on the chosen mapping between 
off-shell and on-shell phase spaces. However, the induced numerical differences remain within 
the accuracy of the calculation, i.e. 0{aTw/Mw)- At high energies, say above 2TeV, when 
peaks in the cross-section become much more pronounced, a more sizeable implementation 
dependence may occur. The present calculation is primarily meant for LEP2 energies and 
slightly above, say up to 500 GeV. 

On the practical side, the results give an insight in the size of RC for off-shell VT-pair pro- 
duction. Within the claimed accuracy it is a complete calculation. It should be stated that 
imaginary parts of loop diagrams have been neglected in the expectation that they only in- 
duce small effects. Moreover, such terms are characterized by an explicit sin 03 and/or sin 04 
dependence in the cross-section, since they select the antisymmetric parts of the "D-matrices. 
Integration over these azimuthal angles removes the imaginary parts of loop diagrams from 
the cross-section. So they do not contribute to the distributions of Sect. ^. It should also be 
stated that some large corrections (ISR, FSR), which usually require the inclusion of higher- 
order terms, are considered purely in first order here. The corresponding higher-order terms 
can be determined in a straightforward way within our approach, using the usual exponentia- 
tion/factorization techniques. 

We have presented the results for leptonic final states. The reason is that those states are 
theoretically well defined. In exactly the same way also quark final states can be treated as 
long as one assumes certain masses for the quarks. It is this mass assignment which gives 
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some arbitrariness in the actual numbers. The RC presented are corrections to ideal theoretical 
situations, which cannot be realized experimentally in the same way. For that purpose the 
calculation should be implemented in an event generator. In principle this is possible. Events 
can be generated in the on-shell phase space with a radiatively corrected weight. The outside 
Breit-Wigner distributions can then generate the invariant masses and consequently off-shell 
events could be constructed from the on-shell ones (with certain angles kept fixed). Event 
generators offer the possibility to include realistic experimental cuts and therefore to study 
effects like the hne-shape deformation in practical cases. The actual numbers presented here 
should give the reader an indication of the size of RC in ideal cases, of which remnants survive 
in practical situations. 

For some questions the present study could already be useful in its present form. An example 
of this would be the comparison of a DPA Born cross-section with CP-conserving anomalous 
triple gauge-boson couplings and a normal DPA cross-section with RC. The question of how 
RC mimic anomalous couplings could be studied in this way, but goes beyond the size and 
scope of the present paper. 

In conclusion, the DPA method for unstable-particle production has been shown to give 
reasonable results. It could also be applied to other unstable-particle production processes that 
undergo electroweak or QCD radiative corrections. 
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A Helicity amplitudes for the virtual factorizable cor- 
rections 



In this appendix we give the basic ingredients for the calculation of the virtual factorizable 
corrections. The one-loop RC to on-shell W^-pair production have been carried out in the 
literature in terms of helicity amplitudes with a particular choice for the decomposition into 
basic matrix elements and invariant functions |12|. This calculation serves as our basis for 
obtaining the RC to the production density matrix V[x-^x^]ix'_^x'^]{Mw , Mw) , defined in Eq. (p^S]). 
Therefore we will set up our conventions in close analogy to what has been used in the numerical 
routines of Ref. \T2 . Once we have fixed the choice of polarization basis in the production stage, 
the same choice should be applied to the decay stages as well, i.e to Vx^x'XMw) in Eq. (p!6D . 



A.l Virtual corrections to the production stage 



The amplitude Ilxj^x2{Mw, for on-shell VT-pair production depends on the spinors of the 
initial-state and on the polarization vectors (pj, Aj) of the W bosons. In order to define 
the latter we first introduce 



^l,2(Pl,2, II) 



ql^iM^ + qli{M^ + t) + p^^it - u) 



ut - Js - AMI 



i-W 



Aivpr 



ls{ut - M^) 

spf-2M^{q, + q2r 
MwJs{s-AM^) 
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using the conventions defined in Sect. |^ and e'^^^'^ = —1. The helicity states, defined in the LAB 
frame, can be expressed in terms of the energy E and velocity P = p/E = ^1 — / E"^ of the 
W bosons in the following way: 



e1{p„±l) = -j= 



et{pi,0) = etipuL) 
eUp2,±l) 



e^,{p„\\)±^e^i{pi, 
E 



1 



-^(0,-l,±^,0), 



1 

V2 



Mw 

£2{P2,\\)Ti£2{P2 

E 



(A 0,0,1), 



V2 



(0,1, ±2,0), 



(/?, 0,0,-1). 



w 



This polarization basis satisfies the usual identities 



Pi£i,,.{Pt, Ai) = 0, ef{pi, Xi) el^{pi, A-) = - 5a,a;. 



(59) 



(60) 
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In its most general form, the amplitude Il\-^\^{M]v , M\y) can be written as a sum of invariant 
functions Cj{a) multiplied by Lorentz-invariant basic matrix elements A^°'(Ai,A2). The basic 
matrix elements are simple, purely kinematical objects and contain the complete dependence 
on the VT-boson polarizations. The invariant functions contain the dynamical information, 
i.e. details of the model, but are independent of the FT-boson polarizations. Both parts depend 
on the helicity of the electron, Ag- = a/2 {a = ±1). In view of our massless treatment of the 
e^, the helicity of the positron is fixed to Ae+ = — Ag- in the virtual corrections. For a specific 
helicity of the electron the decomposition of the helicity amplitudes reads 

9 

U,-MX,{Mw,Mw) = ^C,((t)A<J(Ai, A2), (61) 

with 

-M?(Ai,A2) 

^2(Al,A2) 
^3(Al,A2) 

-MI(Ai,A2) 

-M?(Ai,A2) 
^^(Ai,A2) 
-M^(Ai,A2) 
The helicity projectors 

UJa = ^ ' ^ '° , (63) 

with 75 = i'y^'-f^'-f^'-f'^ ^ project on right- and left-handed massless fermions. Note that our set 
of 18 basic matrix elements Ai'j is overcomplete. Because of the massless treatment of the 
fermions, CP invariance implies the relation^ 

A1j(Ai,A2) =A<J(-A2,-Ai), (64) 

resulting in only 12 independent matrix elements. The last three pairs, A^ygg, have been kept 
for convenience and can be expressed in terms of the others according to 

M^j = ^aM^s+M^2 + ^-j^M^, M^q = -^MI aM^^ = 2[M1+M1 + Mt) 
(65) 

^We have fixed the overall phase of the matrix elements such that this relation holds. The density matrix 
will of course not be affected by this choice. 



v{qi)T^iUJ„u{q2) iPi£2)iP2£i] 



v{qi 



v{pie2) - h{V2(^\ 



uj^u{q 



2), 



v{qi)h{ii - il)^2^^aU{q2), 



v{q 



^^aUiq2), 



i{qie2) - ^2{q2£ij 
v{qi)^iuj„u{q2) {qie2){q2£i), 
v{qi):^iuj^u{q2)\ [{Pie2)iq2i^i) + (p2£^i)(gi£^2)] , 
v{qi)'y''uj^u{q2)\ e^^pr £2 {Vi - P2Y , 

v{qi)Y^M(l2)]^f,uprP2Pl [^UPi^2) " el{p2ei)]. 



(62) 
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For the lowest-order matrix element only a few of these basic matrix elements are relevant: 
Kmx2 (Mw, Mw) = Ml{\, A2) + 2 G2(s, a) [MI{K A2) - A^?(Ai, A2)] , (66) 

where the coefficients Gi^2 are defined as 

ie^ I ^ , . . ofQe , cw ge{,cr) 



2s^y t 



s sw s — Af| 

Here sw and cw are the sine and cosine of the weak mixing angle 9w- 

Mw 



Cw = COS 9 



w 



Sw = sinO 



w 



1 - r2 



and geicr) denotes the coupling between the Z boson and electrons: 



^7e( + l) = -— Qe, ge{-l) 
Cw 

where Qe = ~1 is the charge of the electron. 



1 + IQf.Sy^ 

2swCw 



(67) 



(6^ 



(69) 



For the virtual corrections to on-shell VT-pair production we need the complete list of basic 
matrix elements ■M'j- For the TV-boson helicity states (|59|) the non-vanishing matrix elements 
read: 



A^i(±l,±l 
MI{±1,±1 

A^fi(±l,±l 



MI{±1,0) 




MI{±1,0) 


= Ml{0,Tl) 




= MUO,Tl) 


-^6 (±1,0) 


= MUO,Tl) 


A<?(±1,0) 


= MUO,Tl) 




= -M8(0,Tl) 



2E^{coseT(^) sine 
-2E^{coseT(r) sine 
E^(3 sin^ 9 

2E^f3 sine 
2E'^ {cos e - 13) sine 
-2E^ cos ^ sin ^ 
E^(3 sin^ e 

V2E^. 



(70) 



(71) 



Mw 
y/2E 
Mw 
V2E 

Mw 



2E'/3{coseTcr) 



E' 



2/3 - 2 cos ^ T a{l - p^) (cos e t (r) 



E^{p + 2 cos e±cr) (cos ^ ^ ct) 



E'^p (P + cose) sin^e 



w 

y^2E^P^sinH 
Mw 

^2E'p\cose: 
Mw 
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(0,0) 




2E^P{1 + P^)sme 


A^2(0,0) 


E^ 


8E^/?3 sin e 


-^3(0,0) 


E^ 


8E^f3sme 






[3/?-/53-2cos^] sin^ 


A^5(0,0) 




4:E^{P + cos 9) sine 


-^6(0,0) 




2E'^(3 (/3 + cos^)2sin^ 


A1^(0,0) 


E^ 


8E^(3'^{f3 + cose)sme. 



(72) 



(73) 



From this list and the invariant functions of Ref. |]T2| the density matrix V[\-^^\2][\[x'^]{Mw , Mw] 
follows in a straightforward way. 



A. 2 Virtual corrections to the decay stages 

Since we have chosen a specific polarization basis for the calculation of the production stage, 
both at lowest order and at virtual one-loop order, the same basis has to be used for describing 
the on-shell IV-boson decays. All results presented in this subsection are therefore given in the 
LAB frame, rather than the often used rest frame of the decaying W boson. 

Like in the on-shell VT-pair-production case, we again write the decay matrix element as a 
sum of invariant functions £'j^'* multiplied by Lorentz-invariant basic matrix elements Ai^^\Xi): 

^'C\Mw) = E^rX-''^(Ai), ^t\Mw) = E^i~X"^(A2). (74) 
j j 

In the most general case of the decay oi a W boson into massive quarks, there are four basic 
matrix elements [|15|. For the decay of the W~ boson, W~{p2) — > /2(^2)/2(^2)) they are given 

by 

Mi-\X2) = [u{k2)u;^v{k',)]ie*2k2), 

Mi-\\2) = [uik2)uj+vik',)]ie;k2), (75) 
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with similar expressions for the decay of the boson. For massless particles in the final state 
only Mi"^^ occurs. At lowest-order we then obtain 

Al{Mw) = '-^Mi^\x.), (76) 

with Vfi f the quark mixing matrix. 

For the decay density matrix Vx^x'_{My[r) it is useful to have the expressions for 

M^^\\,K)= E Ml^\X.)M^^^*{K). (77) 

fcrmion helicities 

For the decay of the W~ boson one finds 



sin^ ^4 
(1-/5 cos ^4)^' 



^1 -/3 cos ^4)2^ 



. _§e--^,^;^. ,78) 

The expressions for the charge-conjugate process, describing the decay of the boson, can 
be obtained through the simple relation 

M''o^\Xi,X[) = Mio\-Xi,-X[), with 04 ^ 03, cos 514 cos 6*3, sin 614 -sin 6*3. (79) 

These expressions can be combined with the invariant functions from Ref. |T^ to obtain the 
decay density matrices including virtual RC. It can be seen from Eq. (|78| ) that the matrices 
M-l^^Xi, X'j) contain asymmetric imaginary parts proportional to sin 03 4. These terms will be 
responsible for picking up imaginary loop effects present in the invariant functions, which do 
not depend on 03 4. The symmetric parts of M.f^\Xi, X'^) are real and depend on cos 03^4. Upon 
integration over the azimuthal angles 03,4 the matrices A4^(^\Xi, become real and diagonal, 
and the same holds for the corresponding decay density matrices. 



B The Weyl — van der Waerden formalism for the real- 
photon factorizable corrections 

Like in the case of the virtual factorizable corrections, also for the real-photon factorizable cor- 
rections our choice for the polarization basis and the calculational scheme is guided by existing 
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calculations in the literature. As mentioned in Sect. there is no objection against having 
different choices for the polarization basis in different contributions to the RC, provided that 
the contribution to the density matrix is calculated consistently within the chosen approach. 
We adopt the conventions of Ref. and calculate the real-photon RC in the Weyl-van der 
Waerden formalism. 



B.l The Weyl — van der Waerden formalism for massive gauge bosons 

Before giving the results for the various matrix elements, we first give a few essential details 
of the Weyl-van der Waerden formalism for massive gauge bosons. We follow the conventions 
of Ref. and define the two-dimensional Weyl spinor for a massless particle with light-like 



momentum q as 



_ /ETL^ . 'u- M'. (80) 



{qi + n2)l\/qo -gs 

The indices can be raised and lowered by the spinor metric 



according to 

The spinor products 



^ ~ 1 _1 / ~ -^AB 



q^ = qBe'''', q^ = cab q"" ■ (82) 



{qk) = qAk"" = - q'^kA, {qk)* = q^k"" = - q'^kj, (83) 

are hence antisymmetric. In the Weyl representation for the 7-matrices we obtain the following 
set of translation rules into two-dimensional spinor language: 



M(g,+) 


= v{q,-) 


u{q,-) 


= v{q,+) 






Q^'^Ab 


= Qab 




= qAB = 



q 





A 



<A 







2QK = g^^i^AiJ, 



q^qs 2qk = \{qk)\\ (84) 

with Q, K arbitrary Lorentz vectors and g, k light-like ones. The Dirac spinors u(g, ±) denote 
right-handed (+) and left-handed (— ) states. The matrices a^^^ = ((t°,(t) consist of the 2x2 
unit matrix and the standard Pauli matrices a* [i = 1, 2, 3). 
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For a photon with momentum k we use the polarization vectors describing the two hehcity 
eigenstates 

^ " A_B/_i\^ /9___ 

{kb) ' S I ^^^^^ 

To handle the massive bosons we first decompose their massive momenta pi^2 into a sum 
of two light-like momenta: 



= -r(-l) = v^?i4r=.r(+l). (85) 



Pl,2 = P3,4 + Cl,2 "^3,4, Ci (2^37713) = C2 (2^47714) = M^. (86) 

Note that the so-defined light-like vectors 7713^4 can be chosen freely. An orthogonal basis for 
the three physical polarizations of the massive bosons is now given by 

A B A B 1 

c^Bf,^\ _ , /77 Pa "^3 ^ABf n_,/9 "VPs c-^Bfr,\_^_f^ r, ^ 

4^,^!) _ y§i!±!!5l, (-1) = v^i!!±PL, ,iB(o) = J- (p. - , (87) 

(P4"^4) (^4^4) 

with 1) = £•^'^(+1). It should be stressed that the so-obtained polarization basis does 

not correspond to the helicity eigenstates. However, the corresponding states transform like 
helicity eigenstates under a parity transformation, which is very useful for practical calculations. 



B.1.1 An example: lowest-order on-shell W-pair production 

As an example we apply the above method to the lowest-order on-shell production stage. To 
this end we choose m3 4 = ^4 3 in Eqs. (|86|) and (pT]), and write Ci = C2 = M^/(2p3P4) = c. 

The complete Born amplitude of the process is of the form 

K-Mx^i^w. Mw) = Gi{t) Mtia, X,, A2) + ^2(5, (t)Ms{(t, Ai, A2), (88) 
where the coefficients Gi 2 are defined in Eq. (|67|), and 

Mt{(7,\i,\2) = v{qi)^i{^2-h)h^-u{(l2), 

Ms{(y, \i,\2) = v{qi)f{pi,ei,p2,e2)ujau{q2), (89) 

where 

V''ipueuP2,e2) = 2et{e2Pi) - 2e^^{e^p2) + {p'i - pt){eie2). (90) 

These matrix elements can be translated into two-dimensional representation, e.g. for cr = — 1 
one obtains 

A^t((T = -1, Ai, A2) = gi^^iBAfe -P2)^%iic7 9?' 

;W,(a = -1, Ai, A2) = VBA(Pi'^i'P2,e2)gf • (91) 
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For a — +1 we can make use of the relations 

Mt(a = +1, Ai, A2) = 0, Ms(c7 = +1, Ai, A2) ^ M^a ^ -1, -Ai, -A2), (92) 

where the last identity is the result of parity conservation of the s-channel matrix element 
Ms, since all parity violation is contained in the coefficient G2{s, a). Prom CP invariance one 
obtains two more relations: 

Ms,t(Qi^ -^^ Q2, Ai,p2, A2) = - Mlt(q2, -a, qi, a,p2, -X2,Pi, -Ai), (93) 
so only 6 independent polarization states remain. The independent matrix elements read: 



M(-l,+l,+l) 

M(-i,+i,-i) 
M(-i,-i,+i) 
M(-i,+i,o) 

M(-i,-i,o) 
M(-i,o,o) 



{(llPi){(l2Pi)* 2?iP3 - 



P3P4 

(eiP3)^(?lP4)*(?2P4)* 



P3P4 



V^{qm){q2P3)* 

Mw{P3Pi) 

V2{qiP3){q2PAy 



Mw{pzPaY 



[2q2P4 - 2cq2P3 + 
[2q2P4 - 2cq2P3 - 



(giP3)(?2P3)*| 

for the i-channel matrix elements, and 

Ms{-l,+l,+l 
Ms{- 

Msi-l,-l,0) 

M,(-1,0,0) 
for the s-channel matrix elements. 



1-C + 



:i + c) 



92^4 - cq2Pi 



1,+1,-1) = Ms{-1,-1,+1) = 0, 
V2 



c), 



Mw 



w 



(9lP4)(g2P3)*(P3P4)*(l - C^), 
(9lP3)(P3P4)(g2P4)*(l - C^), 



(?lP3)(g2P3)*Q + 3 - 3C - C^^ 



(94) 



(95) 



B.2 Non-collinear photon radiation from the production stage 

Using the above example as guideline, we now address the process of non-collinear real-photon 
radiation from the production stage: 

e+(gi, ^i)e-{q2: ^2) ^ W+{pu X,)W~{p2, X2HK A). (96) 
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Since we are dealing with non-colhnear radiation and massless initial-state electrons and posi- 
trons, we can ignore the possibility of helicity flip in the initial state. Therefore again the 
condition ai = —02 = —o apphes. As a first step we extend the list of kinematical invariants 
of Sect. 13: 



s = (gi + ^2) , t = {qi -pi) , t 
s' = {P1+P2Y, t' = (g2-P2)^ 
The complete matrix element can now be written in the form 



Wl -P2) 
(^2 - Pi) 



U^iMw, Mw) = -e\G,it')Mj + Gi(t)A<] + ^2(5', a)Ml + ^2(5, a)M 



(97) 



(9^ 



where the functions Gi 2 are the same as the ones defined in Eq. (^). The basic matrix elements 
A^J are invariant under gauge transformations of the radiated photon. They are given by 



M'( = v{q,)\^ 



iizl^ (J. 



2qik 



Ml = v{qi 



2)>t-2 



-2^i{e^Pi) + 2^^{eik) - 2lt{e^ei) 



2pik 



■^2 \ ^^-u{q2), 



+ 



2q2k ' ''^'2p2k 
^7%7-rf(^i'^i'^'2,£2) - Y{pi,ei,P2,e2 



Ml 



viqi)' 



v{qi 



2^2(£^7P2) - H^[e2k) + 2^{e^e2)\ | u;-M(g2) 



^^{6162) + fb{Pl,£l,P2,£2) + fb{P2,B2,Pl,£l)\uJaU{q2 



(99) 



where we introduced the shorthand notation 



Pik 



Vl'{pi,ei,p2,e2] 



+ 4 



- ^ V^'ipi + k, euP2, 82) + [p^ieie^) + e^;{e,k)] (eaPi + e2k) 
sis^){kpi + kp2) - {eik){e^p2) + p^ (£2^7) (^lA;) + {£i£i){,e2Pi) ■ (100) 



The vertex function V can be taken from Eq. (pO|) . Note that the term —2}/^{e^ei) between the 
square brackets of M\ originally had the form (^1 — ]^){e^ei). The difference (^1 + ^)(e^ei) 
cancels against similar terms in Ml- This cancellation is a consequence of the lowest-order 
Ward identity of the boson. In the same way also the Ward identity of the W~ boson has 
been used to simplify Ml- 



For the calculation in the Weyl-van der Waerden formalism we choose 7723 4 = A; in Eqs. (|86D 
and (|R7|). Furthermore we choose the free gauge parameter h in Eq. ( |S5| ) to be equal to qi. Like 
in the case without photon radiation, we can exploit some symmetry relations. First of all CP 
invariance implies the relations 

-^1(91, 92, Ai,p2, A2, fc, A) = Ml*{q2, -cr, gi, a,p2, -A2,Pi, -Ai, k, -A), 

■M3,4(9i,-cr,?2,fT,pi,Ai,p2,A2,A;, A) = A^^^4(g2, -0^, gi, o-,p2, -A2,Pi, -Ai, A;, -A). (101) 

The matrix elements for right-handed electrons are again completely determined: 



Ml{a = +l) = Ml{a 
MUa = +l,K\2,\) = M 



7*1 
3,4' 



a 



-- 0, 

-Ai, — A2, —A), 



(102) 
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where the last identity is the resuh of parity conservation of the s-channel matrix elements 4. 
Due to the symmetry (antisymmetry) property of the quartic (triple) gauge boson vertex under 
the exchange of the W'^ and W~ bosons, one can derive two more relations: 

Ml{pi,\i,p2,X2,k,\) = - A<^(j92, A2,pi, Ai,/c, A), 
M2{pi,XuP2,X2,k,X) = +M2{p2,X2,Pi,Xi,k,X)- (103) 

After all these preparations we now list the independent matrix elements for o" = — 1. In 

order to keep the results as compact as possible we use the shorthand notations {riVj) = {ij) and 
(rjrj) = (ij), with = (91, ?2,P3,P4, ^) for i — (1,2,3,4,5). For the amplitude A1i'(Ai, A2, A) 
we find: 

(15) (24)^ 



MJ{+1,+1,+1) 
MU+l,-l,+l) 
MJ{-1,+1,+1) 
A17(-1,-1,+1) 

M7(0,+l,+l) 

MT(o,-i,+i) 
M7(o,o,+i) 



= -4V2 



(35)2(45) 



(35)2(45)* 
^ (15)(45) 



(35) - (13) - (15) 

(24) ((13) (23)* + (15)(25)*) - C2(13)(45)(35)* 
. (15)" 



(35) 



= -x/2 



(13)2(24)(25)*2 
(15)(45)*(35) ' 



[(13) (23)* + (15) (25)*] [t' + 4(24)] 



(35)2MvK 

2c2(13)(45)(24)*(35)*^, 

(13)2(25)* 
(15)(35)MvK 
(13) (24)* 



(35)(45)Mvt. 

(13)(25)* 
{35){A5)*Mw 
V2{13) 



+ 



t' + 4(24) 



2c. 



t' + 4(24)J, 

{-(13) + (l-c0[(35)-(15)]}, 

^^^^ ((24)(23)* - C2(45)(35)*) + (1 - ci)(24)(25)^ 
(45) (24)* (35)* 



(15) 
(13) 
(15) 



/(13) 
1(15) 



(23)* + (1 -ci)(25)* 



MH+l +1 -1) - ,/o (15)-(13)*-(24)* 

(15)(13)*(25)* 



A17(+1,-1,-1) 
A<I(-1,+1,-1) 
A^7(-l,-l,-l) 



~ (15)*(45)*(35) 

8^2(24)* 
~ (45)(15)*(35)*2 

4v^(24)(25)*2 
(35)*2(45)*(15)^ 



(24)(23)*-C2(45)(35) 
(15) - (35)] [(35) - (13) - (15) 
(35) - (13) - (15) 
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M7(+i,o,-i) 

MJ{0,+1,-1) 
A47(0,0,-l) 



^|||^{-2c.(45)(24)*(35)* + (^^^ 

ri' + 4(24)][(35)-(13)-(15)], 



(15)*(35)*2Mw^ 

4(15)(13)*(24)* 
~ (45)(15)*(35)*Mw. 
2(25)* 



Ci[(15)-(35)] -(13)-(15) + (35)|, 

(15)*(35)*(45)W.l^^^^^^^2^^^^^^*^2^^* 
+ 2((24)(23)* - C2(45)(35)*) [(35) - (13) - (15) 
V2 



(15)*(35)*M2, 



-4c2(24)*(35)*(45) (35) - (13) - (15) 



+ 



i' + 4(24) (ci(15)(13)*(25)* + 2(23)* (35) - (13) - (15) 



The independent matrix elements A4l{\i, A2, A) read: 

4v^(15)(34)' 



M?(+l,+l,+l) 
MI(+1,-1,+1) 



(25) (35) (45) 
2V2 



(25) (35) (45)* [ (15) 



(12) + (35) + (45)], 

(12) 

2)^(45) (14) (23)* + C2(12)(45)(25)*(35)^ 



- (35)^ 



(34) ((12) (23)* + (15) (35)*) - 2(45) (14) 



A^I(-1,-1,+1) 
Ml{0,0,+1) 



2V2(12)=^(34)(25)*^ 
(15)(25)(35)*(45)* ' 

~^ ^ -C2(12)(45)(25)*(34)^ 



(25)(35)Mvi. 



+ 



(104) 



(12)(23)* + (15)(35)*] [(34) + ci(45) - C2(35)] j, 

4(12)(25)* r_c^(i5)(34)(45)* + (13) [(34) + ^^(45) - C2(35)] |, 



(15)(25)(35)*MvK 

-2||(14)(34)*[(34) + C2(35) + ci(45)] 



a/2 r .(13) 



(25) Mf 



w 



+ (13)(35)*[2c^(35) - 4c2(45) + (2c2 - 4ci)(34) - (3ci - 2c2)M^ 
- (14) (45)* [2c? (45) - 4c^(35) + (2ci - 4c2)(34) - (3c2 - 2ci)M^]|. 



(105) 



For Mj{Xi, A2, A) we obtain: 
A^l(+1,+1,+1) = -4%/2 



(15)^^(12)* 
(35)2(45)2 



(34) + (35) + (45) 
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M2(+i,-i,+i) 

A^I(-l,-l,+l) 
A^l(+1,0,+1) 



Ml{-1,0,+1) 



A^I(0,0,+1) 



^^|^{c,(15)(34)(25)*(35)* + 2M(34)(23)*(35) 

+ 2(14) (25)* [(34) + (45) 

V2 (12)(34)^(25)*^ 
2 (35) (45) ' 

4 f(15) 
Mvi.(35)2\(45) 



(34)(23)* (34) + (35)(1 - Cs) + (45)(1 + Ci) 



- (15)(25)*[(34)(1 + 2c2) + (35)(1 + Ca) + (45)(1 + Ci) + M, 

— x |-C2(14) (35) (34) (45)* 
Mvy(45)(35)l ^ ^ 

+ (13)(34) [(34) + (35)(1 - cs) + (45)(1 + ci + 2c2) 
2^2(34) (15) (25)^ 



M2,(35)(45) 



-2ci(l + Ci)(45) + C2(l + C2)(35) 



+ (c2-2ci)(34)+(c2-^ci-^)m, 



2^2(14) 
(45) 



(25)* + (24) 



*(34) 
(35) 



(34) + (35)(l + C2) + (45)(l + ci; 



(106) 



B.2.1 Lowest-order decay of the W bosons 

Having fixed the polarization choice for the real-photon factorizable corrections to the produc- 
tion stage, we now calculate the lowest-order decay parts accordingly, since they are needed for 
obtaining the DPA limit of the full matrix element Aio in Eq. (^). The matrix elements for 
the l^-boson decays are given by 



\/2sw 



Using again 7713,4 = k m. Eqs. (|8^) and (^71) , one ends up with 



Mt\+l) = V2 



{psk)* (psk) 

Mi;\o) = -^{{k[p3){k,p,y - c,{k[k){k,ky 



for the W'^ boson, and 



M}y\+i) = V2 



{k2PA){k'^ky 



{Piky 



V2 



{k2k){k'^p^y 
{P^k) 



(107) 



(108) 
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Mi-\o) = -Lnk,p,){k',p,y - C2{hk){k',ky) (109) 

for the boson. 



B.3 Non-collinear photon radiation from the decay stages 



Next we address the process of non-colhnear real-photon radiation from the decay stages. We 
start off with the decay of the boson: 



iV-ipuXi)-- Mk^)f[{k'Mk,\). 



;iio) 



We do not exphcitly write the hehcities of the final-state fermions. The final-state fermions 
are treated as being massless, hence for non-collinear radiation their helicities are fixed by the 
left-handed interaction with the W bosons: A/^ = — A/^ = +1. The matrix element for process 
(|110|) can be written as 



w 



with 



7(+) 



uik[) 



u{k[) 



a/2 s VI/ 
2kik 



7{+) 



;iii) 



.v{ki 



.,AJ ih + k) 



BC 



2k,k ^-^Dch, 

ik[ + kf^ 



D 



k'.^'e 



-yBA 



2 A^]^ kt 



'IDC 



D 



u{ki) -— r '-LJ-V{ki, 

2pik 



k'i 



lA 



2pik 



;ii2) 



Here the vertex function V can be taken from Eq. (|90|) , but {f, +^1) can be replaced by as 
a result of the lowest-order Ward identity of the boson. 

For the calculation in the Weyl- van der Waerden formalism we choose the same polarization 
basis as adopted for the on-shell W-pair example in Sect. p.l.l|, i.e. 4 = P43. For the 



definition of the photon polarizations we choose the free gauge parameter h in Eq. (|85|) to 
be equal to fci. A straightforward calculation gives the following results for the amplitudes 
-M?^+^(Ai,A): 



-2 



{k[p^){p,ky 



{kik){p3P4)*' 
{k[pA){P3ky 



-2 



{kik){p:iPi) 
V2 



{kik)Mw 



c{k[p,){p,ky-{k[p3){p3ky 



(113) 
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The corresponding expressions for Mli^\Xi,X) read: 





:+!,+!) = 












-1,-1) = 




'^(0,+l) = 




"^(0,-1) = 



-2 
2 
2 



(kik) {k[k) (psp^)* 
c{pm){hPA)*'^ 



{kik)*{k{k)*{p3PAy 
{kik[){k[p,){k^psr 
{k,k){k[k){p;p,) ' 



{kik)*{k[k)*' 
V2{k,k[) 



Mw{kik){k[k) 



{k[p3) (kips)* - c {k[pi) {kiPiY 



{P3P4){kiP3)*{km)* 
Mw{kik)*{k[k)* ' 



(114) 



and 



2pikMl'^^\+l,+l) 
2p,k Ml^^\+l, -1) 
2pikMl^^\-l,+l) 



2c 



'2{k[p3) 
{kik){p3P4)* 
{k[k){psp4){kip4)*^ 
{kik)* (psp^)* 

2{k[p,) 



{kik[) {kip4)* {k[k)* — 2kik{p4ky 



2pikM 



7(+)/ 



-h-V. 



{kik){p3P4) 

ik[k){k,p,Y^ 



{k^k[){kmT{k[kY -2kMP3ky 



2pikMl^^\'d,+l 



2pikMl^^\<d,-l) = 2V2c 



{kik) 
V2 



{k,k[){k[ky{{k[ps) {k^p^r - c {k[p4) {k^p^y 



{kik)Mw 
2k^k[{k[p3){p,ky - c{k[p4){p4k) 
{k[k){p3P4){kip3y{kiPiy 



Mw{kiky 



(115) 



The expressions for the charge-conjugate process, describing the decay of the W boson, can 
be obtained as follows: 



A(-)(M, 



w , 



V2s 



w 



Qf,M 



(116) 



where 



M]^-\X2,X)=[M]''^\-X2,-X)] , with {k,,k[,p3,P4)--{k2,k'„p4,P3). (117) 

When the above matrix elements for real-photon radiation from the decay stages are com- 
bined with the lowest-order matrix element for the production stage, presented in App. [B.l.l| , 
one obtains the DPA limit of the full matrix elements Ai± in Eqs. ( ^6]) and (0). 
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B.4 Radiation of collinear photons 

Up to now we have only discussed the case of non-colhnear photon radiation, which allowed us 
to neglect the fermion masses and the possibility of spin- flip in the initial state. The picture 
changes, however, if the radiated photons are sufficiently collinear with one of the external 
fermions. In such cases factorization takes place, i.e. the matrix element squared including 
collinear radiation can be approximately written in terms of the lowest-order matrix element 
squared and collinear factors. 

Let us first consider collinear photon radiation in the direction of one of the light fermions in 
the production stage of the process, e.g. the positron. In that case the matrix element squared 
can be written in the following form 0: 

J2 l-^coii,e+(gi,(^i,g2,(^2,fc, A)P ^ f!:2\qi,(Ti,k) \M^pj^{xiqi,-a2,q2,(r2)\'^, (118) 

A 

where xi = {E — ko) /E is the ratio of the positron energy after and before photon radiation, 
(71,2 are the helicities of the e^, and 

The last term in this collinear factor gives rise to the so-called spin-fiip, which allows the positron 
to have the same helicity as the electron. Note that we have only indicated the momenta and 
helicities of the relevant particles (e^,7) and that the photon helicities are summed over, as 
the photon cannot be detected anyway. Collinear radiation in the direction of the initial-state 
electron can be obtained in the same way, with the role of the e"*" and e~ interchanged. If the 
initial-state particles are not polarized, as is the case at LEP2, the collinear factor takes on the 
well known form 

When the photon angles are integrated out, the terms oc 1/qik yield contributions of the large- 
logarithmic type [oc log(s/mg)], whereas the term oc ml/{qik)'^ gives rise to additional 0{1) 
contributions, which would have been neglected in a massless treatment of the initial state. 

In the case of collinear photon radiation in the direction of one of the final-state fermions, 
say the fermion /2 from the W~ decay, the factorization reads 

l-^coii,/2(fci, k[, fca, k'2, k)\^ ^ e^Q). 

where the summation is performed over all final-state helicities and = E^^/ {E4 + /cq) is the 
ratio of the /2 energy after and before photon radiation. The other final-state collinear factors 
can be obtained in the same way. 



1 + xj 1 1 + xf ml 
Xi{l — Xi) qik 2xi {qiky 



+ 5, 



(o"liO'2) 



[1 - Xi)^ ml 
2xi {qiky 



(119) 



l+yl 1 



m. 



y2 k2k {k2ky 



\M'j,PAikuk[,k2/y2,k',)\\ (121; 
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C Special integrals for semi-soft photon radiation 



In this appendix we have a closer look at the inclusive treatment of the photon in shifted Breit- 
Wigner resonances l/[-Di + in the vicinity of the Mf resonance (see Sect. ^). We start off 
with factorizable real-photon radiation, involving the ratios 



I A 



I A + 2kpi 



1 



1 



D* + 2kpi Di + 2kpi 



\D, 



IDA 



2iMwTw 



Re 



Di + 2kpi 



(122) 



In order to study the phenomenon of hard-photon suppression we consider the generic integrals 



\Di 



Re 



dk 



l\rn-2 



\s<ko<A 



(27r)32A;o kS A + 2kp^ 



(123) 



for n = 1 or 2. The integration is performed over the photon angles and the photon-energy 
range Xg < ko < A, where is a soft-photon cut-off (A^ -C Tw)- For n = 2 this integral 
quantifies the influence of the shifted resonance on the Mf distribution in the vicinity of the 
pole = M^. For n = 1 it quantifies the effect of 0{k) shifts in the definition of the DPA 
residues. In the latter case we find 



\D, 



Im < Li2 



-1+/3 



Li-; 



-1-/5 



(124) 



with z = Di/{2EA) and Li2 the usual dilogarithm. One can immediately read off that Ji is 
suppressed by 0(Tw/E), irrespective of the precise value for A. For n = 2 the integral reads 



Im 



D* 



MwTw 



1 + log 



/ Di 



\2EX. 



+ 7777— \og{z+l-p) — — \og{z+l + P) 



2P 



2P 



(125) 

This type of integral will lead to an 0{1) contribution. The dependence on the cut-off A, 
however, is suppressed by 0(r^/A^). So, the more energetic the photon is the more suppressed 
its contribution will be. Hence, as soon as A is taken to be much larger than F^y it can safely 
be replaced by infinity. 

Based on the latter observation, we can now list the relevant integrals needed for the inclusive 
treatment of final-state radiation effects involving shifted Breit-Wigner resonances (see columns 
3,4 of Table |T] and column 3 of Table 0). For the radiation from the "''-boson decay stage the 
following four integrals are required: 



dk 



1 



1 



ko>Xa 



fco>As 



{2n)^2ko {2kkif[Di + 2kpi] IGn^mj^Di 

dk 1 _ 

(27r)32A;o i2kki){2kpi)[Di + 2kpi] " 



log 



1 



DiE^ 



+ Li2 1 



1 



log 



D1E3 
+ U2 



log 



. m 



h 

1±£E_' 
2 Es, 
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dk 



1 



1 



fco>As 



ko>Xs 



dk 1 
(27r)32A;o {2kki){2kk[)[Di + 2kpi 



log 



( 2XsMw 



1 



(27r)32fco (2fcpi)[Di + 2A;pi 



1 



1-/3 ' 
1 



(126) 



From these integrals one can determine the correction factor corresponding to the |X_^ | term in 

Eq. my- 



dk 



ko>Xs 



(27r)32A;o 

+ Q5,|iog 
+ Q?j(iog 



^Re 

IT 



log 
log 



ml 

ml, 
■ii 



log 



2XsMw_ 



'-2j3''^[TTp 



+ Li2 1 



1 -/3 E 



2 



+ Li2 1 



Li2 1 



1-P E 



2 E' 



^ +Li2 1 



1+13 E 
2 Es 

1+P E 
2 E', 



(127) 



Here E'^ denotes the energy of f[, i.e. E'^ = E — E^. The correction factor corresponding to the 
|Xi| term in Eq. (H) is obtained by replacing (Di, /i, /(, E^, E'^) by (D2, /2, /2, -^4, -E4). 

Finally we study the hard-photon suppression for the non-factorizable corrections. To this 
end we consider the integrals 



J„ = Re 



dk 



(27r)32A;o k^ A + '^kpi 



:i28) 



^As<fco<A 

for n = 1 and 2. From the results for /i^2 one straightforwardly obtains 



WE(3Mw 



Re <^ 



Lig 



-1 + /3 



Li2 



-1-/? 



(129) 



and 

J2-- 



47r2 



Re n + log 



2E\. 



z+l-p 
2(3 



\og{z + l-l3) 



z + l + (3 
2(3 



\og{z + l+(3)\. (130) 



Again a suppression of 0{Ty// E) is observed for = 1, whereas for n = 2 the dependence on 
the cut-off A is suppressed by 0{Mw'^w I [E . So, again A can be replaced by infinity if it 
is sufficiently large. For explicit expressions for the non-factorizable corrections we refer to the 
literature % 
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